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Abstract 

Aspects  of  a  constitutive  model  for  characterizing  crystalline  metals  containing  a  distribu¬ 
tion  of  dislocation  and  disclination  defects  are  presented.  Kinematics,  balance  laws,  and 
general  kinetic  relations  are  developed  -  from  the  perspective  of  multiscale  volume  averaging 
-  upon  examination  of  a  deforming  crystalline  element  containing  a  distribution  of  displace¬ 
ment  discontinuities  in  the  form  of  translational  and  rotational  lattice  defects,  i.e.,  dislocations 
and  disclinations.  The  macroscopic  kinematic  description  is  characterized  by  a  three-term  mul¬ 
tiplicative  decomposition  of  the  deformation  gradient.  The  micro-level  description  follows 
from  an  additive  decomposition  of  an  affine  connection  into  contributions  from  populations 
of  dislocations  and  disclinations  to  the  distortion  of  the  lattice  directors.  Standard  balance 
equations  apply  at  the  macroscopic  scale,  while  momentum  balances  reminiscent  of  those 
encountered  in  micropolar  elasticity  (i.e.,  couple  stress  theory)  are  imposed  at  the  micro-level 
on  first  and  second  order  moment  stresses  associated  with  geometrically  necessary  defects. 
Thermodynamic  restrictions  are  presented,  and  general  features  of  kinetic  relations  are  postu¬ 
lated  for  time  rates  of  inelastic  deformations  and  internal  variables.  Micropolar  rotations  are 
incorporated  to  capture  physics  that  geometrically  necessary  dislocations  stemming  from  first 
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order  gradients  of  elastic  or  plastic  parts  of  the  total  deformation  gradient  may  alone  be 
unable  to  reflect,  including  evolution  of  defect  substructure  at  multiple  length  scales  and 
incompatible  lattice  misorientation  gradients  arising  in  ductile  single  crystals  subjected  to 
nominally  homogeneous  deformation. 

©  2005  Elsevier  Ltd.  All  rights  reserved. 

Keywords:  Constitutive  behavior;  Crystal  plasticity;  Disclinations;  Dislocations;  Microstructures 


1.  Introduction 

Processes  of  grain  subdivision  and  dislocation  substructure  formation  are  known 
to  substantially  affect  slip  system  activity,  strain  hardening,  stored  lattice  energy,  and 
texture  evolution  in  single  crystals  and  polycrystals  (Hughes  et  al.,  1997,  2003;  Kuhl- 
mann-Wilsdorf,  1999;  Butler  et  al.,  2000;  Barton  and  Dawson,  2001;  Hughes,  2001; 
Leffers,  2001).  Also  measured  within  pure  ductile  metals  and  certain  alloys  at  large 
deformations  and/or  high  temperatures  are  long  range  internal  stress  fields  associ¬ 
ated  with  misoriented  subgrain  boundaries  (Gibeling  and  Nix,  1980;  Argon  and 
Takeuchi,  1981);  moreover,  recent  experiments  on  cyclic  dislocation  cell  structures 
(Kassner  et  al.,  2002)  suggest  such  internal  stress  fields  may  be  of  negligible  magni¬ 
tude  in  the  macroscopically  unloaded  configuration  in  FCC  crystals.  The  trend  in 
many  metallic  systems  of  increasing  strength  with  decreasing  size  of  considered  vol¬ 
ume  or  microstructural  features  is  also  well  known.  Perhaps  most  often  cited  in  this 
context  is  the  Hall-Petch  relation,  in  which  yield  stress  and  cleavage  strength  in¬ 
crease  with  decreasing  grain  size  in  polycrystals  (Hall,  1951;  Petch,  1953).  Since  poly¬ 
crystalline  grain  boundaries  and  deformation-induced  internal  boundaries  share 
certain  fundamental  characteristics  such  as  lattice  incompatibility  (i.e.,  misorienta¬ 
tion  across  interfaces)  and  nominally  large  defect  densities  relative  to  the  bulk  mate¬ 
rial,  one  would  expect  some  degree  of  similarity  in  the  effects  engendered  by  these 
boundaries  upon  homogenized  mechanical  properties  of  the  material  such  as  macro¬ 
scopic  strength  and  hardness. 

Composite  models  wherein  walled  cellular  dislocation  structures,  represented  by 
hard  and  soft  regions  of  relatively  high  and  low  defect  densities,  respectively,  have 
been  suggested  to  explain  the  effects  of  evolving  populations  of  these  substructures 
on  flow  stress  (Mughrabi,  1983,  1988,  2001;  Berveiller  et  al.,  1993;  Zaiser,  1998). 
Meyers  and  co-workers  (Meyers  and  Ashworth,  1982;  Benson  et  al.,  2001;  Fu 
et  al.,  2001)  invoked  composite  models  featuring  grain  boundary  layers  of  relatively 
high  dislocation  density  to  explain  grain  size  influences  on  yielding.  Formulations 
have  also  been  proposed  that  explicitly  embed  subdivision  and  related  dislocation 
substructure  effects  into  the  kinematics  of  crystal  plasticity  theory  (Leffers,  1994; 
Butler  and  McDowell,  1998)  and  the  hardening  and  intergranular  interaction  laws 
of  polycrystalline  models  (Horstemeyer  and  McDowell,  1998;  Horstemeyer  et  al., 
1999).  Composite  micromechanics  models  have  also  been  applied  to  describe  the 
breakdown  of  the  Hall-Petch  effect  in  nanocrystalline  solids  (Capolungo  et  al., 
2005). 
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A  great  many  continuum  treatments  of  elasticity  or  elastoplasticity  include 
enhancements  intended  to  describe  influences  of  microstructure  on  the  mechanical 
response  of  the  material,  perhaps  beginning  with  the  pioneering  director  theory  of 
Cosserat  and  Cosserat  (1909).  Multi-polar  and  higher  gradient  elasticity  theories  fea¬ 
turing  higher  order  stresses  (e.g.,  couple  stresses)  were  developed  by  Truesdell  and 
Toupin  (1960),  Toupin  (1962,  1964),  Green  and  Rivlin  (1964),  Teodosiu  (1967), 
and  Eringen  (1972),  among  others.  Kroner  (1963)  explained  the  origin  of  couple 
stresses  from  a  volume  averaging  perspective,  with  the  Cauchy  stress  for  a  represen¬ 
tative  volume  containing  uniformly  distributed  lattice  defects  defined  as  the  average 
of  the  microscopic  stress,  and  the  couple  stresses  produced  by  fluctuations  in  the 
microscopic  stress  fields  induced  by  the  defects.  Early  strain  gradient  plasticity  the¬ 
ories  were  postulated  by  Fox  (1968),  Teodosiu  (1969),  Lardner  (1969),  Dillon  and 
Kratochvil  (1970),  and  Dillon  and  Perzyna  (1972).  Holt  (1970)  and  Bammann  and 
Aifantis  (1982)  developed  gradient-based  models  of  dislocation  kinetics.  Kroner 
(1973)  and  Hartley  (1975)  considered  the  effects  of  higher  order  moments  of  disloca¬ 
tion  distributions.  By  statistical  averaging  of  elastic  stress  fields  imparted  by  distri¬ 
butions  of  discrete  dislocation  dipoles,  Groma  et  al.  (2003)  derived  continuum 
expressions  for  resistance  to  plastic  deformation,  with  flow  stresses  (due  to  shorter- 
range  dislocation  force  interactions)  depending  upon  the  average  dislocation  density 
and  back  stresses  (due  to  longer-range  interactions)  depending  upon  the  first  spatial 
gradient  of  the  dislocation  density.  In  recent  years,  higher  order  plasticity  theories 
have  exhibited  a  resurgence  in  popularity,  as  evidenced  by  an  abundance  of  models 
featured  in  the  literature  (cf.  Aifantis,  1987;  Fleck  and  Hutchinson,  1993;  Nagdhi 
and  Srinivasa,  1993;  Fleck  et  al.,  1994;  Le  and  Stumpf,  1996a;  Shizawa  and  Zbib, 
1999;  Acharya  and  Bassani,  2000;  Acharya,  2003;  Bammann,  2001;  Regueiro 
et  al.,  2002;  Gurtin,  2002;  Svendsen,  2002;  Garikipati,  2003;  Clayton  et  al.,  2004a). 
The  fundamental  kinematic  descriptions  invoked  by  many  of  these  theories  are  based 
upon  differential-geometric  foundations  often  credited  to  Bilby  et  al.  (1955),  Gunther 
(1958),  Kroner  (1960),  Kondo  (1964),  Truesdell  and  Noll  (1965),  and  Noll  (1967), 
while  the  concept  of  compatibility-maintaining  geometrically  necessary  dislocations 
was  established  by  Nye’s  tensorial  description  (1953)  and  then  popularized  by  Ash¬ 
by’s  (1970)  seminal  paper.  Higher  order  plasticity  theories  have  been  used  to  explain 
grain  size  effects  upon  hardness  or  flow  stress  (Shu  and  Fleck,  1999;  Evers  et  al., 

2002) ,  to  predict  geometrical  characteristics  of  deformation  bands  (Aifantis,  1987; 
Gurtin,  2000),  to  provide  insight  into  dislocation  nucleation  criteria  (Acharya, 

2003) ,  and  to  alleviate  mesh  dependency  issues  in  numerical  implementations 
(cf.  Menzel  and  Steinmann,  2000). 

It  should  be  noted  that  in  the  majority  of  the  above  gradient-type  theories,  geo¬ 
metrically  necessary  dislocations  representing  defect  substructure  will  not  evolve  in 
single  crystals  subjected  to  nominally  homogeneous  boundary  and  initial  conditions; 
in  these  theories,  stress  and  elastic  deformation  gradient  fields  will  remain  spatially 
constant  within  regions  of  a  single  crystal  subjected  to  such  boundary  conditions  un¬ 
less  the  constitutive  formulation  favors  development  of  deformation  heterogeneity, 
for  example  from  the  standpoint  of  material  instability.  By  homogeneous  boundary 
conditions  we  mean  linear  displacement  or  constant  traction  applied  to  the  external 
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surface  of  the  single  crystal,  along  with  any  higher  order  boundary  conditions  -  often 
necessitated  by  the  presence  of  strain  gradients  in  the  formulation,  for  example  -  that 
would  result  in  null  lattice  curvature  on  the  boundary.  Correspondingly,  homoge¬ 
neous  initial  conditions  imply  no  initial  lattice  curvature  or  gradients  in  initial  hard¬ 
ness  within  the  crystal.  Kratochvil  and  Orlova  (1990)  and  Sedlacek  et  al.  (2001) 
addressed  self-organization  of  dislocation  substructure  using  nonlinear  stability 
analyses.  Ortiz  and  co-workers  (Ortiz  and  Repetto,  1999;  Ortiz  et  al.,  2000)  con¬ 
structed  a  nonlocal  lamination  theory  for  modeling  additional  internal  degrees  of 
freedom  associated  with  dislocation  substructure,  with  kinetics  dictated  by  minimi¬ 
zation  of  the  incremental  work  of  deformation  arising  from  a  non-convex  strain  en¬ 
ergy  potential  that  admits  multiple  quasi-stable  states.  Similarly,  Carstensen  et  al. 
(2002)  remarked  how  heterogeneous  plastic  flow  may  result  from  a  lack  of  convexity 
of  the  free  energy  when  evolution  of  deformation  is  viewed  from  the  standpoint  of  a 
constrained  energy  minimization  problem.  Possible  sources  of  instability  or  non¬ 
convexity  of  the  strain  energy  density  of  single  crystals,  from  the  perspective  of  con¬ 
tinuum  plasticity  theory,  include  geometrical  softening  and  strong  latent  hardening 
(Ortiz  and  Repetto,  1999;  Sedlacek  et  al.,  2001;  Carstensen  et  al.,  2002).  Such  lack  of 
convexity  is  thought  to  promote  subdivision  of  crystals  into  subdomains  character¬ 
ized  by  single  slip  and  varying  gradients  of  deformation,  analogous  to  the  minimum 
energy  configurations  postulated  for  microstructural  refinement  in  the  crystallo¬ 
graphic  theory  of  martensite  (Ball  and  James,  1987;  Bhattacharya,  1991). 

In  the  present  work  we  appeal  to  the  mathematical  concept  of  a  disclination  to  de¬ 
scribe  rotational  defects  in  crystals.  In  a  general  sense,  such  defects  may  include  low 
angle,  high  angle,  and  twin  boundaries.  In  essence,  disclinations  embody  equivalent 
arrays  of  dislocations  with  net  Burgers  vectors  that  can  also  be  used  to  represent 
these  boundaries.  Limiting  these  boundaries  to  finite  strain  energy  necessitates  that 
they  be  composed  of  disclinations  of  finite  size,  i.e.,  partial  disclinations  bounded  by 
walls  of  disclination  dipoles,  termed  disclination  structural  units  by  Nazarov  et  al. 
(2000).  The  principal  reason  we  seek  to  extend  finite  elastoplasticity  to  incorporate 
disclinations  is  that  such  defects  convey  important  implications  with  regard  to  size 
effects  and  self-organization  of  dislocation  substructure  which  are  becoming  increas¬ 
ingly  evident  in  experimental  characterization  at  finer  length  scales  (Pantleon,  1996; 
Hughes  et  al.,  1998;  Kassner  and  Perez-Prado,  2000;  Valiev  et  al.,  2002).  By  includ¬ 
ing  both  dislocations  and  disclinations  in  our  framework,  we  are  able  to  distinguish 
between  geometrically  necessary  defect  densities  reflecting  lattice  incompatibilities  at 
multiple  length  scales  (i.e.,  multiple  characteristic  spacings).  For  example,  Hughes 
and  co-workers  (Hughes  et  al.,  1997,  1998;  Hughes  and  Hansen,  2000;  Hughes, 
2001)  have  observed  within  deforming  ductile  FCC  metals  the  formation  of  cells 
of  relatively  small  misorientation  organized  collectively  into  larger  cell  blocks,  with 
average  misorientations  between  blocks  usually  significantly  greater  in  magnitude 
than  those  between  cells.  Upon  increasing  applied  strain,  cell  block  sizes  generally 
decrease  at  faster  rate  than  do  cell  sizes  (Hughes  et  al.,  1997).  In  the  context  of 
our  theory,  the  disclination  concept  can  be  used  to  capture  the  gradients  of  lattice 
rotation  at  the  cell  block  boundaries  that  arise  from  the  organization  and  super¬ 
position  of  relatively  small  misorientations  between  the  cells,  reflected  here  by 
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geometrically  necessary  dislocations.  Additionally,  when  the  kinetics  of  evolution  of 
statistically  stored  defects,  geometrically  necessary  dislocations,  and  geometrically 
necessary  disclinations  are  properly  coupled,  cells  and  cell  blocks  will  emerge  in  sin¬ 
gle  crystals  upon  homogeneous  loading,  as  observed  in  the  aforementioned  experi¬ 
ments,  and  the  subdivided  crystal  will  attain  an  energetically  favorable 
configuration  (i.e.,  a  local  minimum  in  free  energy  over  its  entire  volume).  We  sug¬ 
gest  that  a  lack  of  local  convexity  -  or,  more  precisely,  lack  of  “cross-quasiconvex- 
ivity”  in  the  terminology  of  Carstensen  et  al.  (2002)  -  stems  from  the  superposition 
of  free  energy  wells  associated  with  different  mechanisms,  in  our  case  associated  with 
generation  and  interaction  of  defect  densities  of  various  origins  (e.g.,  populations  of 
geometrically  necessary  and  statistically  stored  dislocations  and  disclinations).  This 
concept  is  demonstrated  explicitly  in  Appendix  A  of  the  present  paper. 

Some  brief  descriptive  remarks  regarding  disclination  defects  are  now  in  order.  In 
the  early  twentieth  century,  Volterra  (1907)  introduced  six  fundamental  types  of  de¬ 
fects  in  elastic  bodies:  three  types  of  translational  displacement  discontinuities, 
known  as  edge  and  screw  dislocations,  and  three  types  of  rotational  incompatibili¬ 
ties,  later  termed  disclinations  by  Frank  (1958)  and  further  classified  as  either  wedge 
or  twist  disclinations  (see  Fig.  1).  Disclination  theory  has  been  applied  to  numerous 
problems  of  interest.  These  include  descriptions  of  micropolar  rotations  in  liquid 
crystals  (cf.  Frank,  1958;  Cermelli  and  Fried,  2002),  rotational  defect  substructures 
and  commensurate  strain  hardening  in  metal  forming  processes  (Romanov,  1993; 
Seefeldt  and  Klimanek,  1997,  1998;  Valiev  et  al.,  2002),  grain  boundary  structure 
in  Bravais  crystals  (Li,  1972;  Gertsman  et  al.,  1989;  Nazarov  and  Romanov,  1989; 
Rybin  et  al.,  1993),  deformation  twins  (Milliner  and  Romanov,  1994),  and  polycrys¬ 
talline  triple  junctions  (Bollmann,  1991).  Disclinations  have  also  been  recognized  as 
characteristic  defects  in  polymers  (Li  and  Gilman,  1970)  and  nanocrystals  (Nazarov 


Fig.  1.  Volterra’s  defects:  (a)  reference  cylinder  with  defect  line  |0  and  cut  surface  S';  (b,c)  edge  dislocations 
and  (d)  screw  dislocation  with  Burgers  vector  b;  (e.f)  twist  disclinations  and  (g)  wedge  disclination  with 
Frank  vector  m. 
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et  al.,  1993;  Konstantinidis  and  Aifantis,  1998).  Molecular  dynamics  simulations 
incorporating  disclination  concepts  (Shenderova  and  Brenner,  1999;  Nazarov 
et  al.,  2000)  have  been  undertaken  to  characterize  grain  boundary  energy  distribu¬ 
tions  over  a  range  of  intergranular  misorientations. 

Continuum  theories  of  distributed  dislocations  and  disclinations  can  be  found  in 
the  geometrically  oriented  papers  of  Anthony  (1969),  Eringen  and  Claus  (1970), 
Lardner  (1973),  Kossecka  and  de  Wit  (1977),  Minagawa  (1977,  1979,  1981),  Amari 
(1981),  and  De  Wit  (1981).  The  texts  of  Nabarro  (1967),  Lardner  (1974),  Mura 
(1982),  Maugin  (1993),  and  Zubov  (1997)  also  include  precise  mathematical  descrip¬ 
tions  and/or  elastic  solutions.  De  Wit  (1973)  developed  a  theory  of  distributed  discli¬ 
nation  loops.  Pcchcrski  (1983,  1985)  employed  disclination  concepts  in  continuum 
formulations  describing  finite  elastoplastic  kinematics,  strain  hardening,  dislocation 
substructure  development,  and  geometrical  softening,  the  latter  resulting  from  local 
lattice  rotations  and  acting  as  a  potential  precursor  to  shear  localization.  An  exten¬ 
sive  review  of  disclination  theory  focusing  upon  defect  kinetics  and  contributions  to 
plastic  strain  hardening  was  provided  by  Seefeldt  (2001),  who  suggested  that  partial 
disclination  dipoles  be  used  to  describe  the  collective  effects  of  dislocations  compris¬ 
ing  misoriented  subgranular  interfaces  (i.e.,  cell  block  boundaries),  manifesting  dis¬ 
location  substructure  refinement  in  advanced  stages  of  plastic  deformation.  Seefeldt 
et  al.  (2001a, b)  used  dislocation-disclination  models  to  predict  texture  diffusion  and 
cellular  refinement  commensurate  with  grain  subdivision.  Panin  (1998)  and  Makarov 
et  al.  (1999)  emphasized  the  role  of  disclination  structures  accompanying  finite 
inelastic  rotations  at  the  mesoscale,  often  occurring  in  conjunction  with  shear  local¬ 
ization  triggered  by  microscopic  heterogeneity.  Lazar  and  Maugin  (2004)  recently 
described  the  stress  field  of  a  wedge  disclination  via  higher  gradient  elasticity  theory. 

Classification  of  disclinations  as  fundamental  defects  in  Bravais  crystals  has  been 
an  occasional  subject  of  debate  in  the  literature  (Kroner,  1983;  Marcinkowski,  1990; 
Kroner  and  Lagoudas,  1992).  This  argument  arises  because  one  generally  must  con¬ 
sider  a  larger  length  scale  of  observation  (e.g.,  the  vector  r  from  the  axis  of  rotation 
to  the  point  of  discontinuity)  to  define  the  net  displacement  jump  associated  with  a 
single  disclination  in  a  Bravais  crystal  (Kondo,  1964;  Romanov,  1993)  than  is  re¬ 
quired  for  an  isolated  dislocation,  whose  effects  manifest  at  a  scale  on  the  order  of 
the  lattice  parameter.  In  fact,  one  may  generally  construct  an  equivalent  lattice  con¬ 
figuration  by  replacing  disclinations  with  organized  arrays  of  translational  disloca¬ 
tions  (Li,  1972;  De  Wit,  1981;  Kroner,  1983;  Marcinkowski,  1990;  Seefeldt,  2001). 
While  such  a  reconstruction  in  terms  of  dislocations  may  be  capable  of  capturing  ex¬ 
actly  the  stress  fields  associated  with  the  original  distribution  of  disclinations,  the  lat¬ 
tice  curvature  may  not  always  be  so  well  reproduced  (De  Wit,  1981).  It  should  be 
noted  that  partial  disclinations  are  needed  to  describe  lattice  misorientations  of 
strength  less  than  90°  in  a  cubic  lattice  (60°  in  a  hexagonal  lattice),  and  disclination 
dipole  descriptions  are  generally  required  to  allow  termination  of  misoriented  inter¬ 
faces  over  finite  distances,  providing  for  bounded  elastic  energies  (De  Wit,  1971, 
1981;  Seefeldt,  2001). 

Pig.  2  demonstrates  how  a  simple  tilt  boundary  may  be  represented  either  in  terms 
of  dislocations  or  partial  disclination  dipoles  (Li,  1972).  In  Pig.  2,  tp  is  the  angle  of 
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Fig.  2.  Tilt  boundary  (left)  comprised  of  edge  dislocations  (center)  or  partial  disclination  dipoles  (right). 


misorientation,  b  is  the  magnitude  of  the  Burgers  vector,  co  is  the  magnitude  of  the 
Frank  vector  (i.e.,  the  strength  of  the  individual  disclination),  h  is  the  spacing  be¬ 
tween  edge  dislocations,  /  is  the  spacing  between  disclination  dipoles,  and  2 r  is  the 
characteristic  spacing  between  partial  wedge  disclinations  comprising  each  disclina¬ 
tion  dipole.  Pure  twist  boundaries  may  be  constructed  from  either  screw  dislocations 
or  partial  twist  disclination  dipoles  in  an  analogous  fashion  (Seefeldt,  2001),  and 
general  grain  boundaries  -  including  all  types  of  coincident-site  interfaces  such  as 
deformation  twins  -  can  be  built  up  from  a  mixture  of  dislocation  or  disclination 
types.  A  certain  amount  of  stored  energy  is  associated  with  the  stacking  fault  pro¬ 
duced  by  each  partial  disclination,  analogous  to  the  stacking  fault  energy  associated 
with  partial  dislocations.  Li  (1960,  1972)  demonstrated  for  an  isotropic  linear-elastic 
solid  the  equivalence  between  the  strain  energy  density  of  a  wedge  disclination  dipole 
and  a  finite  wall  of  edge  dislocations.  We  note  that  for  general  boundaries,  while  a 
suitable  representation  of  the  lattice  configuration  at  the  microscopic  level  in  terms 
of  discrete  dislocations  alone  is  always  possible,  we  find  the  disclination  concept  a 
useful  vehicle  for  engendering  additional  degrees-of-freedom  to  our  continuum  the¬ 
ory  to  distinguish  and  isolate  the  role  of  long  range,  cooperative  defects  in  accommo¬ 
dation  of  lattice  strain  and  curvature. 

In  the  present  paper  we  develop  a  continuum  theory  of  finite  elastoplasticity  for 
crystalline  metals  containing  continuous  distributions  of  dislocation  and  disclination 
defects.  Our  theory  distinguishes,  in  a  novel  manner,  between  kinematic  and  ener¬ 
getic  contributions  from  geometrically  necessary  dislocations  attributed  to  heteroge¬ 
neous  plastic  deformation  fields  and  disclination-type  boundary  structures 
associated  with  initial  and  evolving  subgranular  lattice  orientation  gradients.  A 
three-term  multiplicative  decomposition  of  the  deformation  gradient  is  suggested 
in  Section  2  following  examination  of  a  deforming  crystalline  volume  element 
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containing  ensembles  of  lattice  defects.  A  linear  connection  describing  spatial  varia¬ 
tions  in  the  configuration  of  a  triad  of  lattice  director  vectors  is  introduced,  consisting 
of  contributions  of  stretch  and  rotation  from  the  macroscopic  lattice  deformation 
field  and  micropolar  degrees-of-freedom  that  capture  additional  lattice  rotation  due 
to  disclinations.  Standard  methods  of  differential  geometry  are  employed  to  construct 
the  density  tensors  of  geometrically  necessary  dislocations  and  disclinations,  with  the 
former  accommodating  incompatible  elastic  (or  plastic)  deformation  gradients  asso¬ 
ciated  with  heterogeneity  of  dislocation  glide,  and  the  latter  capturing  additional  mis- 
orientations  across  subgrain  boundaries,  for  example.  In  Section  3,  we  outline 
fundamental  balance  laws,  thermodynamic  restrictions,  energy  density  functions, 
and  kinetic  prescriptions  needed  to  complete  the  theory.  Section  4  consolidates  the 
main  elements  of  the  theory  and  offers  comparisons  with  dislocation-based  plasticity 
models  from  the  literature.  A  discussion  of  potential  applications  follows  in  the  clos¬ 
ing  remarks,  with  a  simple  example  given  in  Appendix  A  illustrating  how  the  frame¬ 
work  can  facilitate  description  of  the  formation  of  cellular  substructures  observed  in 
FCC  metals  at  low  homologous  temperatures  and  large  strains. 

The  following  notation  is  used.  Vector  and  tensor  quantities  are  typically  repre¬ 
sented  with  boldface  type,  while  scalars  and  individual  components  of  vectors  and 
tensors  are  written  in  italics.  The  index  notation  is  often  invoked  for  clarity,  follow¬ 
ing  the  Einstein  summation  convention  and  distinguishing  between  covariant  (sub¬ 
script)  and  contravariant  (superscript)  components,  unless  noted  otherwise. 
Current  configuration  indices  are  written  in  lower  case  Latin,  reference  configuration 
indices  in  upper  case  Latin,  and  intermediate  configuration  indices  are  written  using 
Greek  symbols.  Juxtaposition  implies  summation  over  two  repeated  adjacent  indices 
(e.g.,  (AB)a*  =  AacBch).  The  dot  (scalar)  product  of  vectors  is  represented  by  the 
symbol  (e.g.,  a  •  b  =  a“gabbh,  with  gah  components  of  the  metric  tensor).  Angled 
brackets  denote  a  dual  (scalar)  product  (e.g.,  for  second-rank  tensors, 
<A,B>  =  tr(AB)  =  AahBba,  and  for  contra-covector  pairs,  (a,b)  =  a ab“).  The  colon 
denotes  contraction  over  repeated  pairs  of  indices  (e.g.,  A :B  =  tr(ATB)  =  AahBah 
and  C:A  -  ChcdAcd).  The  symbol  “®”  represents  the  tensor  (outer)  product  (e.g., 
(a®b)o/l  =  a“bh).  Superposed  -1,  T,  and  denote  inverse,  transpose,  and  material 
time  derivative  operations,  respectively.  Additional  notation  is  clarified  later  as  it 
appears  in  the  text. 


2.  Kinematics 

Consider  a  crystalline  volume  element  cut  from  a  single  crystal  (though  perhaps 
originally  embedded  within  a  polycrystal)  and  containing  one  or  more  lattice  discon¬ 
tinuities  (i.e.,  dislocation  and/or  disclination  defects),  as  shown  in  Fig.  3.  We  invoke 
the  following  notation  for  configurations  of  volume  elements:  v0  =  b0cz  B0, 
v  =  bcB,  and  v  =  b  cz  B,  where  global  reference,  intermediate,  and  current  configu¬ 
rations  of  the  entire  macroscopic  body  are  labeled  B0,  B,  and  B,  respectively.  Typical 
dimensions  of  the  referential  volume  v0  are  assumed  to  adhere  to  the  following 
inequalities: 
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(a)  (b)  (c) 

Fig.  3.  Crystal  volume  element  with  defects:  (a)  reference  configuration  b0,  (b)  unloaded  intermediate 
configuration  b,  and  (c)  current  configuration  b. 


«0  rref  ^  Lg, 


(1) 


where  /ref  =  is  the  scalar  effective  length  of  an  edge  of  the  reference  volume  ele¬ 
ment  (Fig.  3(a)),  and  a0  and  La  are  the  lattice  parameter  and  average  grain  diameter, 
respectively.  Taking  a0  «  10~lom  and  Lq  «  10_5m-10_3m  as  representative  of 
engineering  metals,  we  assume  /ref  «  10  3  m  10  4  m  to  satisfy  Eq.  (1).  In  Fig.  3 
the  sizes  of  the  defects  (e.g.,  Burgers  vectors)  relative  to  the  size  of  the  volume  ele¬ 
ment  are  exaggerated  for  clarity  of  presentation.  It  is  noted  that  we  do  not  require 
the  referential  volume  to  constitute  a  statistically  representative  volume  element  in 
the  sense  of  Hill  (1972).  Low  defect  densities  are  admissible.  Dislocations  are  repre¬ 
sented  by  the  symbol  _!_,  and  disclinations  by  the  symbol  >. 

The  three  configurations  of  the  volume  element  shown  in  Fig.  3  are  defined  as  fol¬ 
lows.  The  reference  configuration  volume  v0  is  the  crystal  lattice  as  it  existed  prior  to 
application  of  applied  external  forces  (typically,  and  rather  arbitrarily,  at  an  initial 
time),  such  that  it  is  free  of  traction  along  the  external  boundary  s0.  The  current  con¬ 
figuration  volume  v  is  the  deformed  crystal  lattice,  with  non-vanishing  traction  vec¬ 
tor  t  per  unit  current  area  applied  on  its  external  boundary.  The  defects  present  in  v0 
may  slip  over  a  finite  distance  as  a  result  of  the  applied  loading  and  may  exit  the  vol¬ 
ume  element  in  the  current  configuration,  and  additional  defects  may  enter  the  vol¬ 
ume  element  during  the  history  of  loading.  The  intermediate  configuration  volume  v 
is  achieved  by  unloading  the  current  volume  v  by  an  inverse  elastic  deformation  field 
Fe  “ 1  associated  with  traction  removal  from  the  surface  of  the  volume,  as  will  be  de¬ 
fined  more  precisely  later.  We  remark  that  here  v  is  not  necessarily  the  usual  inter¬ 
mediate  configuration  of  multiplicative  finite  elastoplasticity  (cf.  Lee  and  Liu, 
1967;  Asaro,  1983).  The  local  deformation  gradient  f  for  material  points  with  local 
coordinates  XA  within  the  volume  element  v0  is  defined  as  the  tangent  mapping  be¬ 
tween  reference  and  spatial  configurations,  i.e., 


f  =  T(t>x  =  0  GA  =  x°ASa ' 


(2) 


where  g„  and  G4  are  basis  vectors  and  covectors  assigned  to  the  current  and  refer¬ 
ence  states,  respectively,  and  xa  =  q>a(XA  j)  away  from  displacement  discontinuities. 
We  note  that  f  cannot  be  defined  by  Eq.  (2)  in  locations  in  the  immediate  vicinity  of 
slip  planes  over  which  dislocations  have  passed,  for  example,  where  the  lattice 
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deformation  is  discontinuous.  In  other  words,  it  is  a  discontinuous,  anholonomic 
mapping  in  such  locations  (cf.  Clayton  et  ah,  2004b). 

The  total  deformation  gradient  F  for  the  crystal  element  is  specified  by  the  motion 
of  its  external  boundary,  assuming  homogeneous  basis  vectors  in  the  reference  and 
current  configurations  for  simplicity,  i.e., 


F  = 


1 

Vo 


xanA  dso 


ga  C8>  GA  =  F^a  ®  Gf 


(3) 


with  n  d.s'o  an  oriented  reference  area  element,  as  shown  in  Fig.  3.  When  a  defect  line 
crosses  the  element’s  boundary  in  the  current  state,  the  coordinates  xa  will  be  multi¬ 
valued  functions  of  XA.  In  such  cases  we  take  for  the  contribution  to  (3)  an  average 
of  the  contributions  of  the  coordinates  x“+  and  xa~  on  opposite  faces  of  the  discon¬ 
tinuity  at  the  boundary,  i.e.,  (2v0)  ~1fSo(xa+  +xa~)nA  ds0-  Notice  that  when  the  local 
deformation  held  is  purely  elastic,  meaning  no  defect  generation  or  motion  occurs 
between  reference  and  current  states,  xa  are  continuous,  single-valued,  and  differen¬ 
tiable  functions  of  XA,  and  Eq.  (3)  reduces  to 


K 


l 

^0 


dr0  =  — 
vo 


(4) 


upon  invocation  of  the  generalized  Gauss’s  theorem  (Hill,  1972;  also  see  derivation 
in  Clayton  and  McDowell,  2003). 

We  define  at  the  scale  of  the  volume  element  a  recoverable  elastic  deformation  gra¬ 
dient  Fe  =  VeRe.  The  elastic  stretch  Ve  is  associated  with  the  average  externally  ap¬ 
plied  stress  acting  on  the  element  and  is  determined  explicitly  from  (Clayton  and 
McDowell,  2003) 


Ve  = 


n  ds0 


n  dso 


(5) 


where  x  are  the  local  coordinates  of  the  external  boundary  of  the  element  corre¬ 
sponding  to  a  specially  defined,  second  traction-free  intermediate  configuration  b, 
i.e.,  with  t  =  0,  as  shown  in  Fig.  4.  Configuration  b  arises  from  instantaneous  re¬ 
moval  of  traction  along  the  boundary,  constrained  such  that  rotation  of  the  volume 
element,  Rel,  does  not  occur.  The  unloading  procedure  is  assumed  instantaneous  - 
to  avoid  viscous  inelasticity  upon  load  removal  -  with  inertial  effects  neglected  (Kra- 
tochvil,  1971;  Clayton  and  McDowell,  2003).  Disclinations  and  dislocations  within 
the  volume  element  are  locked  in  during  this  hypothetical  process,  with  no  defect 
motion  or  generation/annihilation  permitted  during  unloading.  Correspondingly, 
the  defect  populations  in  Figs.  3(b)  and  (c)  are  identical  apart  from  the  net  elastic 
unloading  and  rigid  body  motion  of  the  volume  element.  Since  lattice  defects  are 
not  introduced  during  unloading,  Fe  (or  its  inverse)  does  not  alter  the  holonomicity 
(or  lack  thereof)  of  the  material  within  the  volume  element,  at  the  microscale.  How¬ 
ever,  since  the  deformation  achieved  upon  unloading  may  not  be  compatible  across 
neighboring  volume  elements,  Fe  is  generally  anholonomic  (i.e.,  incompatible)  at  the 
macroscale,  in  this  case  defined  as  corresponding  to  length  scales  significantly  larger 
than  that  of  the  volume  element  shown  in  Fig.  3.  In  other  words,  when  traction  is 
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Fig.  4.  Configurations  and  director  vectors. 


removed  from  the  volume  of  the  element,  its  boundary  may  undergo  deformation 
that  is  compatible  or  incompatible  with  neighboring  elements,  with  the  degree  of 
incompatibility  depending  upon  the  scale  of  volume  considered  relative  to  the  distri¬ 
bution  of  defects  contained  within  the  volume.  The  elastic  rotation  tensor  Re  is  deter¬ 
mined  such  that  each  member  of  a  triad  of  lattice  director  vectors  assigned  to  the 
centroid  of  the  crystal  element  deforms  between  configurations  b  and  b  as 

da  =  (F*-iyadx,  (6) 

with  (a,  a  =  1,2,3).  As  depicted  in  Fig.  4,  the  director  vectors  in  the  current  ( b )  and 
intermediate  (. b )  configurations  are  labeled  with  Latin  and  Greek  subscripts  as  d„  and 
d7,  respectively.  Notice  that  these  directors  map  between  configurations  in  an  anal¬ 
ogous  fashion  to  components  of  the  “slip  plane  normal  vectors”  of  conventional 
crystal  plasticity  theory  (Asaro,  1983),  and  that  Re  is  assumed  to  encompass  any  ri¬ 
gid  body  rotations  of  the  entire  (poly)crystal  in  global  configuration  B.  Furthermore, 
when  temperature  changes  occur,  Fe  is  assumed  to  capture  the  net  thermal  expansion 
or  contraction  of  the  volume  element  and  the  corresponding  change  in  the  lattice 
directors,  most  typically  assumed  in  practice  to  be  an  isotropic  deformation  mode 
in  engineering  metals  with  cubic  crystalline  structure,  though  possibly  anisotropic 
in  non-cubic  lattices  (cf.  Nye,  1957).  In  this  context,  the  volume  element  in  thermoe- 
lastically  unloaded  configuration  b  is  of  uniform  reference  temperature  and  null 
average  stress,  and  Fe  is  referred  to  as  the  recoverable  thermoelastic  deformation  gra¬ 
dient.  The  total  lattice  deformation  gradient ,  introduced  as  F^  in  Fig.  4,  is  discussed 
in  detail  later  in  the  text. 

We  next  introduce  the  average  plastic  deformation  gradient  for  the  volume  ele¬ 
ment,  written  here  as  Fp.  As  demonstrated  by  Teodosiu  (1969),  Rice  (1971),  and  Or¬ 
tiz  and  Repetto  (1999),  the  discontinuous  plastic  deformation  gradient  in  the  vicinity 
of  a  single  moving  dislocation  may  be  estimated  in  terms  of  the  Burgers  vector  and 
slip  plane  geometry.  When  multiple  mobile  defects  are  considered,  the  time  rate  of 
plastic  deformation  over  some  finite  volume  (e.g.,  that  of  Fig.  3)  may  be  character¬ 
ized  in  terms  of  the  defect  fluxes,  i.e., 
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Lp  =  FPFp-‘  =  1^2  ®  t  +  E  ®  tj  ■  «>  (7) 

with  Kk  a  scalar  geometry  factor  (Li  and  Gilman,  1970)  and  £  the  rank  three  covar¬ 
iant  permutation  tensor.  Generally  hereafter,  we  use  the  overbar  notation  to  refer  to 
quantities  framed  in  configuration  b  of  Fig.  4.  In  Eq.  (7),  summation  is  applied  over  j 
populations  of  straight  dislocation  lines,  each  with  velocity  vector  |+,  and  k  popula¬ 
tions  of  straight  disclination  lines,  each  with  velocity  vector  £+  and  effective  disclina- 
tion  radius  rk.  The  total  (as  opposed  to  net)  dislocation  and  disclination  density 
tensors  introduced  in  (7)  are  defined  as 

H  =  (K  +  k)K ® ok,  =  +  ii‘K ® t ,  (8) 

with  p{  and  p>_  the  densities  (line  length  per  unit  volume)  of  positively  and  negatively 
signed  dislocations,  where,  for  each  value  of  j,  the  positive  and  negative  dislocations 
share  the  same  tangent  line  but  oppositely  oriented  Burgers  vectors  b/  =  —If  and 
velocities  |+  =  — 1_.  Analogously,  rjk+  and  fjk_  are  densities  (line  length  per  unit  vol¬ 
ume)  of  signed  disclinations  with  tangent  lines  C  and  Frank  vectors  of  =  —mk_  .  Eq. 
(7)  extends  previous  works  of  Teodosiu  (1969),  Lardner  (1974),  and  Werne  and  Kelly 
(1978)  to  include  disclination  line  contributions,  following  the  paper  by  Li  and  Gil¬ 
man  (1970)  on  disclinations  in  polymers.  Similar  expressions  could  be  formulated  to 
include  curved  defect  lines  (cf.  Teodosiu,  1969;  Das  et  al.,  1973),  though  these  tend  to 
be  somewhat  more  complex  as  the  tangent  vectors  of  the  defects  are  no  longer  spa¬ 
tially  constant.  For  example,  Teodosiu  (1969)  derived  an  expression  for  the  velocity 
gradient  contribution  from  a  mobile  density  of  strictly  circular  dislocation  loops.  The 
discrete  approximation  of  (7)  and  (8)  may  be  acceptable  given  that  contemporary  3D 
dislocation  dynamics  models  often  piecewise-linearize  curved  dislocation  segments 
(cf.  Zbib  et  al.,  1998). 

Eqs.  (7)  and  (8)  have  been  introduced  for  illustrative  purposes,  to  show  how  mov¬ 
ing  discrete  crystal  defects  contribute  to  the  time  rate  of  change  of  Fp.  Alternatively, 
we  often  find  it  more  convenient  to  write  Eq.  (7)  from  the  perspective  of  crystal  plas¬ 
ticity  kinematics,  i.e., 

Lp  =  ^fsi0mi  +  Ll,  (9) 

i 

where  summation  runs  over  i  slip  systems,  each  with  shearing  rate  f,  contravariant 
direction  vector  s',  and  plane  normal  co variant  vector  m'.  The  term  Lp  includes  de¬ 
fect  motions  not  readily  describable  in  terms  of  slip  or  pseudo-slip,  for  example  non¬ 
conservative  plastic  flow  stemming  from  dislocation  climb  and  certain  motions  of 
disclination  tangent  lines  and/or  axes  of  rotation  accompanied  by  vacancy  formation 
(Das  et  al.,  1973).  Shearing  rates  on  each  system  may  be  expressed  as 

Y  =  p'bv1  +  VfKrcoz'  +  if  (p  ,  rj'  j ,  (10) 

with  p'  =  p‘+  +  p‘_  the  total  (as  opposed  to  net)  mobile  dislocation  line  density,  b  the 
magnitude  of  the  Burgers  vector,  vl  the  signed  mean  scalar  dislocation  velocity, 
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if  =  Tf+  +  if  the  total  mobile  disclination  density,  co  the  magnitude  of  a  typical 
Frank  vector,  k  and  r  denoting  a  geometric  factor  and  effective  disclination  radius 
both  assumed  identical  over  all  slip  systems  for  simplicity  (this  assumption  can  be 
relaxed  to  address  a  distribution  of  partial  disclination  dipoles),  and  z‘  a  signed  mean 
scalar  velocity  of  disclination  lines  moving  on  slip  system  i,  in  direction  s'.  The  first 
term  in  (10)  follows  from  Orowan  (1940),  the  second  from  Li  and  Gilman  (1970). 
Contributions  of  time  rates  of  generation  and  annihilation  of  mobile  defect  densities 
to  the  slip  rates  are  embodied  by  the  functions  d'  (cf.  Pcchcrski.  1983),  which  may 
additionally  be  extended  to  incorporate  dislocation  flux  through  boundaries  of  finite 
volumes  (including  effects  of  time  rates  of  geometrically  necessary  dislocation  densi¬ 
ties),  as  suggested  by  Arsenlis  et  al.  (2004).  Velocities  are  assumed  positive  for  defects 
moving  in  the  slip  direction  s'  and  negative  for  those  moving  in  the  direction  —s'.  The 
plastic  deformation  gradient  Fp  is  assumed  to  leave  the  lattice  director  vectors  unper¬ 
turbed;  as  shown  in  Fig.  4,  d A  and  d5  are  identical  when  represented  in  a  consistent 
global  coordinate  system.  Notice  also  that  Fp  is  isochoric  when  L(  =  0  and  s'  and  m' 
are  orthogonal,  i.e., 


(11) 


where  Jp  is  the  Jacobian  determinant  of  Fp  and  tr  denotes  the  trace  operation.  We 
remark  that  the  relative  contribution  of  each  of  p‘  and  if  to  the  plastic  velocity  gra¬ 
dient  in  (9)  is  set  forth  implicitly  by  the  choice  of  flow  rule  for  f  (see  later  Eq.  (57)), 
formulation  of  which  may  be  guided  by  heuristic  kinetics  principles  such  as  maximal 
dissipation.  Alternatively,  one  may  approach  the  flow  rule  by  capturing  characteris¬ 
tics  of  dislocation  and  disclination  driving  forces  versus  velocities,  explicitly  describ¬ 
ing  y’  in  (10). 

Consider  for  a  moment  the  usual  multiplicative  decomposition  F  =  FeFp  (Bilby 
et  al.,  1957;  Kroner,  1960;  Lee  and  Liu,  1967)  in  the  context  of  our  volume  averaging 
framework.  We  have  assumed  that  the  entire  deformation  in  fictitious  configuration 
b  is  caused  by  the  time-dependent  generation  and  motion  of  displacement  disconti¬ 
nuities  (dislocations  and  disclinations)  integrated  over  the  activated  slip  planes,  via 
the  definitions  of  Lp  in  Eqs.  (7)  and  (9).  Neglected  in  the  typical  kinematic  descrip¬ 
tion  F  =  FeFp,  when  applied  to  a  crystalline  volume  of  finite  size,  are  the  net  effects  of 
residual  lattice  deformation  fields  contained  within  the  volume  element.  These  local 
residual  deformations  are  attributed  to  the  local  lattice  strain  fields  in  the  vicinity  of 
each  defect  line  and  may  also  include  residual  lattice  deformations  arising  from  ther¬ 
mal  unloading  from  a  heterogeneous  temperature  field.  As  a  generalization  of  the 
two-term  decomposition,  we  propose  a  three-term  multiplicative  decomposition 
for  a  crystal  volume  element  (cf.  Clayton  and  McDowell,  2003;  Clayton  et  al., 
2004a)  consistent  with  Fig.  4,  i.e.. 


F  =  FeF'Fp, 


(12) 


with  F,  Fe,  and  Fp  defined  already,  leaving 
F‘  =  Fe_1FFp_1, 


(13) 
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which  may  further  decomposed  as 

F'  =  ViRi  =  R'U',  (14) 

thus  functioning  as  the  net  contribution  to  F  of  any  residual  elastic  stretch  (V1,  U1) 
and  rotation  (R‘)  of  the  crystal  lattice  remaining  upon  unloading  by  the  applied 
stress.  The  residual  lattice  deformation  arises  from  considering  the  (in)compatibility 
of  microelastic  lattice  strains  and  rotations  within  the  defect  kinematic  fields,  as  per 
the  original  concept  of  Bilby  et  al.  (1955).  The  deformation  (tangent)  map  F1  of  (13) 
is  not  necessarily  a  volume  average,  but  rather  an  indicator  of  the  net  effects  of  local 
deformation  heterogeneity  on  the  total  volume-averaged  deformation  gradient.  Eqs. 
(7)  and  (9)  regard  the  lattice  as  an  aggregate  of  rigid  “blocks”  deforming  by  relative 
sliding  on  slip  planes  between  the  blocks.  The  local  residual  lattice  deformation, 
whose  net  effects  are  embodied  by  F1,  is  necessary  to  deform  each  block  such  that 
the  aggregate  fits  together  into  a  single  continuous  piece  in  configuration  b  (Figs. 
3  and  4),  with  the  only  displacement  discontinuities  within  v  attributed  to  incompat¬ 
ibility  in  the  immediate  vicinity  of  the  defect  lines.  Decomposition  (12)  was  intro¬ 
duced  by  the  present  authors  to  describe  elastoplastic  polycrystals  (Clayton  and 
McDowell,  2003)  and  elastoplastic  single  crystals  (Clayton  et  al.,  2004a).  It  is  intro¬ 
duced  here  from  a  slightly  different  standpoint,  for  (thermo)elastic  single  crystals 
containing  a  distribution  of  dislocation  and  disclination  defects.  We  emphasize  that 
description  (12)  is  scale  dependent,  depending  upon  the  size  of  the  crystal  volume  ele¬ 
ment  to  which  it  is  applied,  and  resolution  dependent,  depending  upon  the  manner  in 
which  local  kinematics  within  the  volume  element  are  resolved  and  accounted  for  in 
the  homogenized  description.  Notice  from  Fig.  4  that  F1  affects  the  representation  of 
the  lattice  directors,  i.e., 

d5=C"da,  (15) 

thereby  accounting  for  slip  system  reorientation  (see  also  later  Eq.  (63))  due  to  de¬ 
fects  contained  within  the  volume  element  in  the  unloaded  state  and  also  any  cumu¬ 
lative  director  motion  due  to  the  disclination  flux  which,  in  contrast  to  the 
dislocation  flux,  generally  induces  a  net  lattice  rotation  (cf.  Fardner,  1973;  Pccherski, 
1983;  and  our  later  Eqs.  (65)-(67)). 

Fig.  5  conceptualizes  the  physics  of  Eq.  (12)  from  the  standpoint  of  (a)  a  Volterra 
process  and  (b)  a  corresponding  lattice  for  a  volume  element  containing  a  single  edge 
dislocation.  Notice  how  the  displacement  discontinuity  in  the  wake  of  the  dislocation 
(i.e.,  left  side  of  the  cylinder  in  Fig.  5(a))  is  sealed  by  the  residual  lattice  deformation 
gradient  F\  Prior  to  the  application  of  F1,  for  the  body  labeled  v  in  Fig.  5,  elastic 
strain  fields  are  absent,  and  the  plastic  deformation  gradient  Fp  is  completely  defined 
in  terms  of  the  relative  motion  of  the  two  halves  of  the  lattice  on  opposite  sides  of  the 
slip  plane  (dotted  line  in  Fig.  5(b)).  Note  also  that  v  is  generally  locally  incompatible 
across  the  entire  slip  path,  containing  slip  discontinuities  but  no  lattice  strains,  and  is 
free  of  internal  residual  stress  fields.  This  is  in  contrast  to  v,  which  may  contain  inter¬ 
nal  residual  stress  and  lattice  strain  fields  but  includes  no  discontinuities  except  for 
those  in  the  immediate  vicinity  of  individual  defect  lines  contained  within  the 
volume. 
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(b) 

Fig.  5.  Configurations  of  crystal  element  containing  edge  dislocation:  (a)  Volterra  model  and  (b)  lattice. 


Please  note  that  the  particular  sequence  of  deformation  maps  in  (12)  is  chosen 
based  on  rational  arguments.  The  plastic  deformation  term  Fp  is  the  rightmost  in 
our  decomposition,  as  it  leaves  the  lattice  vectors  unaltered  (Fig.  4),  permitting  iden¬ 
tification  of  configuration  b  of  Fig.  4  with  Mandel’s  (1973)  isoclinic  configuration. 
The  residual  deformation  due  to  micro-heterogeneity  in  the  presence  of  lattice  de¬ 
fects,  F1,  is  placed  second  in  our  decomposition,  as  it  affects  the  lattice  directors 
via  (15),  yet  is  assumed  unaffected  by  superposed  rigid  body  motion  or  a  change 
in  spatial  coordinate  frame  (Kratochvil,  1972).  Finally,  the  recoverable  elastic  defor¬ 
mation  Fe  logically  assumes  the  leftmost  position,  as  the  stretch  Ve  is  associated  with 
unloading  of  the  average  traction  acting  on  the  volume  element  from  the  current 
state  (Eq.  (5))  and  the  rotation  Re  accounts  for  all  lattice  rotations  not  embodied 
in  F1,  including  net  rigid  body  motions  of  the  solid.  The  rotational  components  R' 
of  F1  and  Rp  of  Fp  are  assumed  to  evolve  independently  of  rigid  body  motions  of 
the  solid.  Thus,  under  a  change  of  spatial  coordinates  x  — >  fix  +  c,  for  which  fi  is 
a  unimodular  rotation  satisfying  Q1  =  il  1  and  c  is  a  spatially  constant  translation 
vector,  the  deformation  maps  transform  according  to  F  — +  OF,  Fe  — +  OFe,  F1  — >  F1, 
and  Fp  ->  Fp. 

We  may  also  choose  to  introduce  a  fictitious  local  natural  configuration  6  free  of 
all  defect  lines,  slipped  regions,  and  commensurate  initial  residual  stresses,  following 
Teodosiu  (1969)  and  Rice  (1971).  Eq.  (12)  is  then  extended  to  jf  =  Fjf0  = 
FcF‘Fp.Jf  o  with  Jf0  :  7T  — >  Tb0  the  deformation  gradient  of  the  volume  element 
due  to  the  creation  of  defect  lines,  grain  boundaries,  etc.  (and  possible  residual  elastic 
lattice  strains)  present  in  the  initial  configuration  v0 ■  For  modeling  the  deformation 
of  a  structure  from  an  engineering  standpoint,  knowledge  of  Jf  0  and  reconstruction 
of  l  are  not  essential  in  our  approach,  since  covariant  components  of  the  total  defor¬ 
mation  gradient  F  are  referred  to  the  actual  state  of  the  body  at  t  =  0,  denoted  here 
by  the  local  volume  element  v0.  Furthermore,  changes  in  free  energy  are  defined  here 
with  respect  to  the  reference  configuration  or  initial  state,  as  opposed  to  the  natural 
state.  The  residual  stresses  and  energy  commensurate  with  Jf0  are  usually  not 
considered  explicitly  in  thermodynamic  treatments  of  finite  crystalline  plasticity, 
even  when  the  map  Jf0  is  introduced  in  the  kinematics  (Teodosiu,  1969;  Rice, 
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1971;  Teodosiu  and  Sidoroff,  1976;  Hartley,  2003).  Also,  we  mention  that  the  afore¬ 
mentioned  theories  refer  the  covariant  leg  of  the  plastic  deformation  gradient  to  the 
generally  anholonomic  natural  state,  arriving  at  a  decomposition  F  =  FeFpJf  0 1 .  This 
contrasts  with  our  model  and  those  of  Le  and  Stumpf  (1996a)  and  Stumpf  and 
Hoppe  (1997)  that  refer  the  covariant  components  of  the  plastic  deformation  gradi¬ 
ent  to  the  (initial)  reference  configuration,  which  is  holonomic  to  the  current  config¬ 
uration.  The  first  fundamental  difference  between  our  stress-free  “natural 
configuration”  6  and  our  stress-free  “intermediate  configuration”  b  is  this:  the  for¬ 
mer  is  achieved  via  unloading  each  crystal  element  from  the  reference  (i.e.,  initial) 
configuration,  while  the  latter  is  realized  via  unloading  each  element  from  the  spatial 
(i.e.,  current)  configuration  (Clayton  et  al.,  2004b).  The  second  difference  is  that 
while  b  is  achieved  by  external  unloading  of  the  local  volume  element,  leaving  this 
volume  intact  at  the  microscale  (as  conceptualized  in  Fig.  4  by  applying  Fe_1  to  b, 
leading  to  configuration  b  with  conditions  t  =  0),  the  natural  configuration  6  can 
only  be  achieved  by  cutting  the  local  volume  element  v0  into  multiple  pieces  such  that 
the  initial  residual  stresses  are  relaxed  at  the  scale  of  the  individual  pieces.  Of  course, 
when  the  initial  configuration  is  a  perfect  lattice,  such  as  in  the  left  side  of  Fig.  5(b), 
free  of  initial  elastic  lattice  deformation  associated  with  defects,  there  is  no  distinc¬ 
tion  between  what  we  call  natural  and  reference  configurations  (i.e.,  6  and  b0  coin¬ 
cide)  and  jf0  =  1.  Although  not  pursued  here,  it  is  possible  to  assign  the  initial 
structure  and  lattice  distortion  fields  corresponding  to  arbitrary  initial  defect  topol¬ 
ogy  as  part  of  .3f  «  if  the  resulting  free  energy  and  residual  strain  fields  are  of  interest. 

Others  have  proposed  three-term  multiplicative  decompositions  for  continuum 
elastoplasticity  exclusive  of  damage,  although  each  within  a  slightly  different  context 
than  the  present  work.  Kratochvil  (1972)  suggested  the  decomposition  F  =  FeF'Fp  for 
polycrystalline  solids,  with  Fe  the  recoverable  elastic  deformation,  Fp  the  viscoplastic 
deformation,  and  F1  accounting  for  any  other  microstructural  rearrangements  -  e.g., 
twinning  and  phase  transitions  -  that  influence  the  physical  properties  and  hence,  the 
specific  free  energy  density.  Horstemeyer  (1995)  and  Horstemeyer  and  McDowell 
(1998)  introduced  the  decomposition  F  =  FeF‘Fp  for  a  single  crystal,  with  Fe  and 
Fp  the  usual  elastic  and  plastic  deformation  gradient  fields  of  crystal  plasticity  the¬ 
ory,  and  with  F1  representing  time-dependent  anelastic  rearrangement  of  the  crystal¬ 
line  lattice,  associated  with  micro-residual  stresses  within  the  crystal,  and  possibly 
accounting  for  the  lattice  deformation  accompanying  grain  subdivision.  Also  consid¬ 
ering  single  crystal  plasticity  theory,  Butler  and  McDowell  (1998)  proposed  the 
decomposition  F  =  FeF'Fp,  with  F1  directly  associated  with  heterogeneous  grain  sub¬ 
division,  distinct  from  the  net  contributions  of  dislocation  glide  embodied  in  Fp.  Lion 
(2000)  introduced  a  decomposition  F  =  FeF’Fp,  with  the  macroscopic  polycrystalline 
non-recoverable  deformation  decomposed  into  residual  elastic  (F1)  and  plastic  (Fp) 
parts.  Bammann  (2001)  and  Regueiro  et  al.  (2002)  suggested  a  similar  three-term 
decomposition  with  F1  the  incompatible  residual  elastic  deformation  gradient  and 
Fe  the  superposed  compatible  elastic  deformation  due  to  the  applied  stress.  Hartley 
(2003)  suggested,  from  a  discrete  dislocation  standpoint,  the  decomposition 
F  =  FeF'FpJf0 1  for  a  single-  or  polycrystalline  volume  element.  In  Hartley’s  model 
(2003),  Fe  is  a  compatible  recoverable  elastic  deformation  gradient  associated  with 
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the  traction  acting  on  a  representative  volume  of  the  material,  and  F1  is  further  mul- 
tiplicatively  decomposed  into  terms  representing  residual  elasticity  due  to  disloca¬ 
tions  contained  within  the  volume  element  and  incompatible  elasticity  due  to 
sources  of  deformation  outside  the  element  (e.g.,  intergranular  incompatibility). 

Our  deformation  mapping  F1  includes  effects  arising  from  defects  contained  within 
the  volume  element  at  time  t,  residual  thermoelastic  strains,  and  internal  boundaries 
and  stacking  faults  left  by  moving  defects  such  as  partial  dislocations  and  partial  dis- 
clination  dipoles.  Lattice-preserving  contributions  from  mobile  defects  that  have  tra¬ 
versed  the  volume  prior  to  time  t  are  embodied  in  the  rate  equations  for  Fp,  (7)  (10). 
It  should  be  noted  that  our  Eq.  (12)  is  a  multiscale  decomposition,  with  the  meaning 
of  each  term  dependent  upon  the  specific  scale  (i.e.,  size  of  (poly)crystalline  volume 
element)  to  which  it  is  applied,  as  discussed  in  Clayton  and  McDowell  (2003).  If  en¬ 
forced  over  a  relatively  homogeneously  deforming  region  encompassing  a  statisti¬ 
cally  large  number  of  grains,  or  a  region  encompassing  a  single  crystal  containing 
a  random  distribution  of  defects  for  which  there  is  no  net  Burgers  vector  (i.e.,  statis¬ 
tically  stored  defects  only),  one  would  expect  F1  — >  1  and  Fe  to  be  a  holonomic  defor¬ 
mation  field  within  this  region,  making  b  a  compatible  configuration.  On  the  other 
hand,  when  (12)  is  applied  to  describe  a  single  crystal  or  polycrystal  subjected  to 
nonuniform  deformation  (i.e.,  strain  gradients),  the  corresponding  volume  element 
will  generally  exhibit  lattice  curvature  due  to  geometrically  necessary  defects,  F1  will 
not  reduce  to  the  identity  mapping  due  to  non-negligible  net  effects  of  residual  lattice 
deformation,  and  Fe  will  be  anholonomic  as  a  result  of  the  heterogeneous  stress  state 
existing  across  neighboring  elements.  We  also  remark  that  continuum  finite  element 
simulations  (Clayton  and  McDowell,  2003;  Clayton  et  ah,  2004a)  suggest  that  the 
magnitude  of  the  stretch  V1  of  Eq.  (14)  remains  small  relative  to  that  of  the  total  ap¬ 
plied  stretch  in  isothermally  deforming  single  crystals  containing  low-angle  subgrain 
boundaries. 

The  deformation  of  the  lattice  director  vectors  between  configurations  b  and  b,  as 
depicted  in  Fig.  4,  is  dictated  by  the  total  lattice  deformation  F:/  =  FeF‘: 

d  a=(F*~l)%  (16) 

Assuming  that  the  directors  in  configuration  b  are  orthonormal,  and  assigning  an 
external  Cartesian  coordinate  system  to  the  anholonomic  space  b  with  corresponding 
metric  tensor 

di  •  d jj  =  gx]j  =  <%;  (17) 

we  define  a  metric  Cy  in  the  spatial  configuration  associated  with  the  total  director 
strain: 


C ^  =  d 

L  ab  ~ 


i  b  = 


_  A  _  Z7i?-la  c  _ 

#  d/i  -  ^  .a  .b  • 


(18) 


Also  introduced  are  coefficients  of  a  linear  (i.e.,  affine)  connection  I  (cf.  Eisen- 
hart,  1926;  Schouten,  1954;  Wang  and  Truesdell,  1973)  describing  the  absolute 
change  of  the  director  vector  field  in  the  current  configuration  in  terms  of  the  covar¬ 
iant  derivative  V,  i.e., 
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Vfedfl  =  -  r£dc,  (19) 

with  the  subscripted  comma  denoting  partial  differentiation  with  respect  to  current 
coordinates  x.  The  connection  coefficients  are  defined  by  (Minagawa,  1979,  1981) 

r-d  -  F%aF%f*  +  Q;l  =  r*  +  q£,  (20) 

where  T'cl  =  F^aFfv'x  =  -F^“cF^~Xci  are  coefficients  of  the  so-called  “crystal  con¬ 
nection”  of  non-Riemannian  dislocation  theories  (Bilby  et  al.,  1955;  Noll,  1967; 
Le  and  Stumpf,  1996b),  and  Q(ah  are  additional  micropolar  degrees-of-freedom  rep¬ 
resenting  contributions  of  distributed  disclinations  to  the  spatial  gradients  of  the  lat¬ 
tice  director  field.  When  disclinations  are  absent  and  the  lattice  directors  are  spatially 
uniform  in  configuration  b ,  the  expression  for  parallel  transport  d„.;,  =  T/:))dr  de¬ 
scribes,  to  first  order,  the  deformed  lattice.  Upon  assuming  that  V*da  =  0  (Minag¬ 
awa,  1981),  the  connection  (20)  allows  one  to  interpolate  for  the  directions  and 
magnitudes  of  the  lattice  directors  between  centroids  of  neighboring  crystal  volume 
elements.  For  example,  in  the  trivial  case  when  =  0,  the  lattice  directors  are  spa¬ 
tially  constant.  The  crystal  connection  (f)  component  of  (20)  accounts  for  effects  of 
first-order  spatial  gradients  of  (the  inverse  of)  V:/ ,  while  the  variable  Q  of  (20)  ac¬ 
counts  for  the  additional  spatial  variations  of  lattice  directors  not  captured  by  the 
first  spatial  gradient  of  F^.  The  average  continuum  deformation  of  the  director  vec¬ 
tors  located  at  the  volume  element’s  centroid  is  determined  by  ,  as  indicated  by 
Eq.  (16).  Lattice  stretch  and  rotation  gradients  (e.g.,  at  the  scale  of  subgrain  cells 
and  cell  blocks)  in  the  current  configuration  are  represented  by  the  coefficients 
(20).  Note  that  unlike  the  deformation  gradient,  I  is  not  a  two-point  tensor  or  con¬ 
figurational  mapping,  but  rather  a  spatial  operator.  The  covariant  components  of  Q 
are  assigned  the  following  anti-symmetry  property  (recall  that  we  use  italic  font  for 
individual  components  of  non-scalar  entities  such  as  Q) 

Qcba  —  Qcb  F da  =  ~Qcab  =  Qc\ba\l  (21) 

with  bracketed  indices  anti-symmetrized,  i.e.,  2gc[fea]  =  Qcba—Qcab ■  Quantities  Qcha 
and  Qcab  are  effectively  equivalent  rotation  measures  only  when  lattice  strains  are  neg¬ 
ligible,  i.e.,  when  Cfa  «  Sda.  From  (21),  the  connection  is  metric,  since  the  covariant 
derivative  of  C :/  vanishes: 


v  cx  _  _  p-“rse  _  p“rse  _0..dcse  _  n..dcse  _  _2n  -  n 

V  c^ab  —  '-'ab.c  1ca'^dh  1  cb^  ad  ^ca^db  SJcb^ad~  A^ic(ab)  ~  u> 


(22) 


=o 


where  parentheses  denote  symmetrization,  e.g.,  2 Qciba)  ~  Qcba  +  Qcab-  Components 
of  the  torsion  tensor  T  of  the  connection  I  are  given  by 

Fcl  =  rci-rb°c  =  T:b  +  2Q^  (23) 

where  T fb  is  the  torsion  of  the  crystal  connection,  often  associated  with  the  density  of 
geometrically  necessary  dislocations  when  disclinations  are  absent  (Bilby  et  al.,  1955; 
Kroner,  1960;  Noll,  1967).  The  components  of  the  Riemann-Christoffel  curvature 
tensor  Rhrad  formed  from  the  connection  coefficients  F£  are  defined  as  (Schouten, 
1954;  Marsden  and  Hughes,  1983) 


228 


J.D.  Clayton  et  al.  /  International  Journal  of  Plasticity  22  (2006)  210-256 


n-.  a  ——  j p..<2  T~<..a  ,  r'..a  r'-.e  r<..a  jp..e 
^ fed  ^  db,c  ^  cb,d  '  ^  ce^  db  ^  de*  cb"> 

which,  because  r£  =  T'£  +  Qcb,  we  are  able  to  rewrite  as 

Kd  =  *£+2V[e0f  +  Q£Q£  -  Q£Qcl  +  T£Q£, 

=0 


(24) 

(25) 


where  the  curvature  from  the  crystal  connection,  R'b'“d,  vanishes  identically  since  f  is 
integrable  (cf.  Le  and  Stumpf,  1996b).  Using  (21),  we  can  express  the  fully  covariant 
rendition  of  the  curvature  tensor,  Rabcd  =  C:/aj  Rh)d,  as  (Minagawa,  1979) 

R[ab)cd  =  2  V[fgdj[Aa]  +  TfdQe [6a],  R(ab)cd  =  2  ^/[cQd\(ab)  +  R cdQe(ab )  =  0-  (26) 

Thus  from  (24)  and  (26),  Rabcd~  R\ah\\,  j\-  and  R  vanishes  completely  when  Q  =  0. 

Consider  a  Burgers  circuit  c  in  the  current  configuration,  enclosing  area  a  com¬ 
prised  of  oriented  differential  elements  n  da.  A  total  Burgers  vector  accounting  for 
the  incompatibility  induced  by  the  torsion  and  curvature  tensors  may  be  written 
as  (Lardner,  1973;  Minagawa,  1979) 

Ba  ee  sdb‘  [  (T£  -  Recy)nd  da  =  B°r  +  (27) 

J  a 

where  2 Bar  =  zdbcf  Tbacnd  da  describes  the  closure  failure  of  c  and  2 B\  =  sbdc 
LRe>e  nd  da  measures  the  change  in  direction  of  position  vector  xe  upon  parallel 
transport  about  c  with  respect  to  the  connection  I We  can  re-write  (27)  in  terms 
of  the  second  rank  geometrically  necessary  dislocation  tensor  «  and  second  rank  geo¬ 
metrically  necessary  disclination  tensor  0,  each  referred  to  the  current  configuration 
b: 


Ba  =  (*“*  +  C*-lafsfgb9gdxb)nd  da,  (28) 

J  a 

where 

2aad  =  g dbcT.M,  4Qgd  =  B*ba  sdce  Rabce .  (29) 

Fig.  6(a)  illustrates  the  total  Burgers  vector  B  =  BT  +  BR  introduced  in  Eq.  (27)  in 
terms  of  parallel  transport  of  a  lattice  director  vector  da  about  an  incompatibility 
circuit  c  (Kondo,  1964).  When  viewed  from  the  standpoint  of  a  Volterra  process 


d  „  a:, 


Fig.  6.  Burgers  vector  with  contributions  from  (a)  dislocations  (BT)  and  disclinations  (BR),  and 
(b)  director  vector  rotation  due  to  60°  wedge  disclination. 
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(Fig.  1),  one  can  imagine  the  body  in  Fig.  6(a)  to  consist  of  the  superposition  of  a 
single  edge  dislocation  (producing  the  incompatibility  BT)  and  a  single  wedge  discli- 
nation  (producing  the  incompatibility  BR).  Fig.  6(b)  shows  the  rotation  of  a  director 
parallel-transported  about  a  complete  60°  wedge  disclination  in  a  hexagonal  lattice 
(Lardner,  1974). 

Notice  that  a  and  9  contain  enough  information  to  fully  reconstruct  T  and  R, 
respectively,  because  of  the  skew-symmetry  properties  of  these  quantities: 

T'cl  =  =  ScbdCCad ,  Rabcd  =  R[ab][cd\  =  ShaeScdf^  ■  (30) 

The  geometrically  necessary  defect  density  tensors  of  (29)  are  related  to  the  summed 
contributions  of  discrete  crystal  defects  as 

a  =  Pi]°j  ®  %  -  6  =  E  ^k(°k  ®  C*>  (31 1 ) 

j  k 

with  pi,  b)  and  S’  the  net  scalar  dislocation  line  density,  Burgers  vector,  and  unit  tan¬ 
gent  line,  respectively,  for  dislocation  population  j,  and  with  rjk,  mk,  and  f*  the  net 
scalar  disclination  line  density,  Frank  vector,  and  unit  tangent  line  vector,  respec¬ 
tively,  for  disclination  population  k.  All  quantities  in  (31)  are  referred  to  the  spatial 
configuration  b. 

Notice  that  a  and  0  do  not  account  for  curved  defect  segments  and  combinations 
of  defect  lines  that  contribute  no  net  Burgers  vector  or  Frank  vector.  This  becomes 
clear  when  we  re-write  (3 1 )  in  terms  of  positively  and  negatively  signed  populations 
(Werne  and  Kelly,  1978): 

a  =  E(^W-)  *4®^  0  =  E('7+-'7-K®Ci,  (32) 

j  k 

where  b(.  =  b  and  cok+  =  - cok_ .  As  discussed  following  (8),  curved  defect  segments 
can  be  treated  via  piecewise  linearization  (Zbib  et  al.,  1998;  Zbib  and  De  La  Rubia, 
2002)  in  order  to  address  the  distribution  of  line  tangents  of  curved  dislocation  lines. 
We  emphasize  that  the  total  mobile  defect  densities  of  (8)  contribute  to  the  plastic 
velocity  gradient  (7),  while  the  net  defect  densities  (both  mobile  and  immobile)  con¬ 
tribute  to  the  equivalent  incompatibility  measures  in  (29),  (31)  and  (32).  Ashby 
(1970)  introduced  the  concept  of  statistically  stored  dislocations  to  represent  defects 
that  accumulate  under  homogeneous  plastic  flow  but  do  not  contribute  to  net  incom¬ 
patibility.  Later  we  appeal  to  this  concept  by  introducing  densities  of  statistically 
stored  dislocations  and  statistically  stored  disclinations  allowing  us  to  quantify  influ¬ 
ences  of  differences  between  total  cumulative  defect  densities  and  the  geometrically 
necessary  densities. 

Linearized  compatibility  equations  for  the  net  defect  density  tensors  of  (29),  (31) 
and  (32)  follow  from  identities  of  Bianchi  and  Schouten  (cf.  Schouten,  1954),  ex¬ 
pressed  in  Cartesian  coordinates  as 

n ^ = R'ihcd]  ->  =  zabc(?c,  R-b;dA  =  o  -  0?  =  o.  (33) 

In  a  small-strain  formulation  (i.e.,  additive  elastoplastic  strains  and  rotation  gradi¬ 
ents),  De  Wit  (1981)  inferred  from  equations  analogous  to  (33)  that  disclinations 
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may  act  as  sources/sinks  for  dislocations,  and  that  disclination  lines  cannot  end 
abruptly  within  a  crystal.  This  linearization,  while  physically  illustrative,  is  not  essen¬ 
tial  in  our  general  finite  strain  and  rotation  framework. 

Complete  pull-backs  of  a  and  0  to  the  intermediate  configuration  b  are  written  as 
follows: 

aals  =  JeFe-l«oiabFeh111 

_  re  rye-lot  bed  ryJ^a  ry£f -la.  rye-lfi  ,  re  rye-la  bde a  rye-l/i  (>XA\ 

—  Jr.a  b  r  .a  r  .d,c  F  .b  ~T~jr.a  b  ^[def  .b  J  \D*J 


~e*p  =  jeF\ 


J  ffb  J7e  1  /' _  je^e-Uadghcen 

o  r  b  ~  &y  C  a  b 


7e-ip 


(35) 


with  f  the  Jacobian  determinant  of  F''.  Variables  a  and  0  may  be  regarded  as  “elas¬ 
tic”  in  the  sense  that  they  are  derived  completely  from  lattice  kinematic  quantities  Fe, 
F1,  Q  and/or  spatial  gradients  of  Fe,  F\  and  Q.  On  the  other  hand,  the  following 
“plastic”  quantities  are  now  introduced  in  the  sense  of  push-forwards  from  the  ref¬ 
erence  configuration: 


^  =  -r'FyscDFyuF^cF';B  - 


-  1A  j ypOt  jyP 
D,Cr  .B 


(36) 


?  -  -rlrA 


=  _  1  y 


<-[GD][CE]f  .B, 


(37) 


with  F  A  =  F‘*FPA  the  residual  deformation  gradient  and  J  the  Jacobian  determinant 
of  F.  By  Q  and  R  we  denote,  respectively,  a  generalized  rank  3  inelastic  micro-rota¬ 
tion  variable  and  rank  4  curvature  tensor  on  the  reference  coordinate  frame  such 
that  the  “total”  defect  densities  aT  and  0T  vanish  in  the  intermediate  configuration, 
i.e., 


aT  =  a  +  otp  =  0,  0T  =  0+0p  =  O.  (38) 

Please  note  that  Q  of  (36)  and  R  of  (37)  are  not  constrained  to  be  compatible  with  a 
linear  connection  defined  on  the  reference  frame,  in  contrast  to  Q  of  (34)  and  R  of 
(35)  which  are  compatible  with  the  connection  coefficients  I '  of  Eq.  (20)  in  the  cur¬ 
rent  configuration.  Equations  given  in  (38)  are  equivalent  to  the  additive  relations  of 
Eringen  and  Claus  (1970)  in  the  linearized  limit.  Also,  when  disclinations  vanish, 
Q  =  0,  0  =  0  =0,  and  the  first  of  (38)  gives  JeFe; 1  1  *Fe~ 1  /f  = 

J  AFPdc^b^  following  automatically  from  relations  F  =  FeF‘Fp  =  F^F11 

and  §c F :/  1  dx  =  §(  Fp  dX  for  current  (reference)  Burgers  circuit  c  (C). 


3.  Balance  laws,  thermodynamics,  and  constitutive  framework 


The  net  nominal  stress  for  the  crystalline  element  illustrated  in  Fig.  3  is  defined  as 


S 


-  /  X  ( 

f’O  Js„ 


1  to  dvo, 


(39) 
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where  t0  is  the  traction  vector  per  unit  reference  area.  Eq.  (39)  reduces  to 
u0S  =  f  s  cit'd,  with  s  the  local  nominal  stress  within  the  volume  element,  only  under 
quasi-static  conditions  and  when  discontinuities  within  the  volume  are  traction-free 
(Nemat-Nasser,  1999).  The  net  Cauchy  stress  is  defined  by 

rh  =J-lFaASAb,  (40) 

with  J  =  dct(F)  \/g l G  the  Jacobian  determinant  of  the  deformation  gradient  F,  and 
g  and  G  determinants  of  metric  tensors  gab  and  GAB  associated,  respectively,  with 
spatial  and  reference  coordinate  systems  shown  in  Fig.  3.  Standard  balances  of  mass, 
linear  momentum,  and  angular  momentum  at  the  level  of  the  volume  element  are 
now,  respectively,  asserted,  i.e., 

Po  =  PJ,  S^  +  pJ*  =  p0xa1  F“ASAb  =  SAaFbA,  (41) 

with  p0,  p,  and  f“  the  reference  mass  density,  current  mass  density,  and  body  force 
component  per  unit  reference  volume,  respectively,  and  the  vertical  bar  denoting 
covariant  differentiation  with  respect  to  the  reference  coordinate  system.  From 
(40)  and  the  last  of  (41),  Sab  -  £(ah\  We  add  the  caveat  that  one  may  choose  to  delve 
at  this  point  into  a  polar  or  hyperstress  theory  (cf.  Toupin,  1962;  Garikipati,  2003) 
for  which  the  Cauchy  stress  tensor  is  non-symmetric  if  dislocation  spacings  are  on 
the  order  of  the  length  scale  of  the  volume  element,  or  if  the  volume  element  has 
a  strong  net  dislocation  density  resulting  in  local  stress  fluctuations  that  may  trans¬ 
late  into  higher  order  moments  of  stress  on  the  same  order  of  magnitude  as  the  aver¬ 
age  stress,  Eq.  (39).  For  simplicity  we  assume  here  that  the  volume  element 
considered  is  large  enough  to  permit  us  to  neglect  higher  order  stresses  in  the  mac¬ 
roscopic  momentum  balances  of  (41),  though  we  shall  account  explicitly  for  the  net 
effects  of  heterogeneity  of  lattice  deformation  through  inclusion  of  F1  and  the  defect 
density  tensors  (representing  lattice  curvature)  in  the  free  energy  function,  micro¬ 
scopic  momentum  balances,  and  kinetics  framework  that  follow. 

Introducing  the  average  internal  energy  per  unit  reference  volume  U0,  effective 
heat  flux  vector  q,  and  average  heat  generation  per  unit  reference  volume  H0,  the  lo¬ 
cal  form  of  the  balance  of  energy  is  written 

U0  =  SAagabFbA-qfA+H0,  (42) 

where  the  quantities  in  (42)  are  averaged  over  a  local  volume  element  of  size  consis¬ 
tent  with  the  preceding  discussion  accompanying  (39).  The  entropy  production  and 
reduced  dissipation  inequalities  are  written  in  the  reference  frame  as 

^  T  “  (  r)  x ’  sAag^A  -  T  ^ T n)  ^  ^  +  W  (43) 

with  N0  the  entropy  per  unit  reference  volume  and  the  Helmholtz  free  energy 
\l/0=  U0—N0T.  Thermodynamic  restrictions  (cf.  Coleman  and  Noll,  1963;  Coleman 
and  Gurtin,  1967)  can  be  deduced  from  inequalities  (43),  the  second  written  in  terms 
of  two-point  tensors  S  and  F  and  the  free  energy  per  unit  reference  volume  We 
find  it  useful  to  express  (43)  in  terms  of  quantities  with  all  components  referred  to 
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intermediate  configuration  b  of  Fig.  3  (i.e.,  the  reference  state  for  instantaneous  ther¬ 
moelastic  response  -  cf.  Scheidler  and  Wright,  2001),  which  we  view  as  the  most  con¬ 
venient  configuration  for  deducing  thermodynamic  restrictions  and  posing 
constitutive  assumptions  (see  also  Regueiro  et  al.,  2002).  After  somewhat  lengthy 
but  standard  manipulations  (see  Clayton  et  al.,  2004a),  the  second  of  (43)  can  be 
written  as 


\  (TO  >^  +  NT,  (44) 

where  Z",  =  is  a  velocity  gradient,  =  JeFe^lLacgcbFeJIb  is  the 

Mandel  stress  (Mandel,  1974),  and  the  heat  flux  vector  referred  to  b  is  defined  by 
q*  =  J  b'‘AqA ,  with  the  anholonomic  temperature  gradient  T  .,  =  T,A  F  r/  .  The  free 
energy  and  entropy  per  unit  volume  in  b  are  found  by  i j/  =J  i//0  and  N  =  J  N 0, 
respectively. 

We  make  the  following  general  assumption  regarding  the  dependency  of  the 
Helmholtz  free  energy  function  for  the  crystalline  volume  element,  referred  to  con¬ 
figuration  b  and  neglecting  the  effects  of  temperature  gradients  within  the  element: 


=  iA(ce,  v',«, 


^5  5  T 1  §  ) 


j. 


(45) 


The  covariant  elastic  strain  tensor  C'),  =  Fe“gabFep  is  included  to  model  the  change  of 
average  elastic  energy  density  with  a  change  of  external  loads,  a  standard  constitu¬ 
tive  assumption  in  finite  crystalline  elastoplasticity  theories  (cf.  Le  et  al.,  1998). 
The  left  stretch  tensor  associated  with  F',  denoted  by  V',  is  incorporated  to  reflect 
contributions  to  the  free  energy  from  residual  microelasticity  within  the  volume  ele¬ 
ment,  and  may  be  non-negligible  when  the  deformation  within  the  volume  element  is 
heterogeneous  (Clayton  and  McDowell,  2003).  The  elastic  energies  due  to  net  lattice 
curvatures  in  the  volume  element  induced  by  geometrically  necessary  dislocations 
and  disclinations  are  reflected,  respectively,  by  the  inclusion  of  a  and  ().  defined  in 
(34)  and  (35).  As  a  point  of  clarification,  recall  that  the  total  mobile  defect  density 
tensors  at)  and  0t  contributing  to  the  plastic  velocity  gradient  were  introduced  in 
(7)  and  (8),  while  the  net  geometrically  necessary  defect  density  tensors  a  and  0  (these 
include  both  mobile  and  immobile  defects)  were  formally  introduced  in  (31)  and  (32). 
Notice  that  when  Q  =  0  (no  disclinations)  and  F‘  -  1,  we  have  Fl!  =  jefeyEabc 
Fe~f  =  Jp-xFpyBCFp(B,  meaning  that  (45)  agrees  with  constitutive  assumptions 
made  in  previous  gradient-based  dislocation  theories  from  the  literature  (cf.  Stein- 
mann,  1996;  Regueiro  et  al.,  2002).  Since,  as  mentioned  already  in  Section  2,  the  ten¬ 
sors  a  and  0  in  (45)  do  not  include  a  measure  of  the  total  length  of  all  dislocation  and 
disclination  lines  within  a  given  volume  element  (examples  include  statistically  stored 
defects  consisting  of  closed  loops  and  tangent  lines  of  opposing  signs),  the  scalar 
parameters  sp  =  b  v/p))  and  i.,t  =  rcoy^  are  incorporated  to  model  the  net  contribu¬ 
tions  of  the  elastic  self-energy  of  the  statistically  stored  dislocation  density  (line 
length  per  unit  volume  in  b),  ps,  and  the  statistically  stored  disclination  density  (also 
a  line  length  per  unit  volume  in  b ),  ijs,  to  the  free  energy.  We  regard  lp  and  as  local 
residual  lattice  strain  measures  associated  with  these  defects  (Bammann,  2001). 
Please  notice  that  heterogeneity  of  lattice  deformation  (reflected  in  (45)  by  V1,  a, 
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and  0 )  is  not  essential  to  engender  substantial  energy  storage,  as  sp  and  ln  account  for 
a  large  fraction  of  the  observed  stored  energy  of  cold  working  under  conditions  of 
homogeneous  plastic  flow  and  increase  in  its  magnitude  with  the  accumulation  of 
the  average  plastic  deformation,  Fp.  Rather,  these  tensorial  variables  tend  to  reflect 
scale  effects,  particularly  the  latter  two.  Additionally,  in  an  effective  sense,  V1  in  (45) 
accounts  for  higher  order  moments  of  distributions  of  lattice  defects  within  the  vol¬ 
ume  element  that  lead  to  internal  heterogeneity,  such  as  dislocation  pile-ups,  for 
example,  that  are  not  accounted  for  by  the  first-order  averages  over  the  volume  ele¬ 
ment  of  defect  fields  a,  0,  sp,  and  er  Reasoning  for  inclusion  of  metric  g  will  become 
clear  later.  Each  variable  in  (45)  is  considered  invariant  with  respect  to  changes  in  the 
spatial  frame  of  reference.  We  choose  not  to  include  the  rotational  part  R1  of  F1  in  the 
free  energy  function  because  material  rotations  of  a  local  volume  element  in  un¬ 
loaded  configuration  b,  while  affecting  orientations  of  the  lattice  directors  via  (15), 
do  not  influence  the  amount  of  stored  elastic  energy  contained  within  that  element. 

We  note  that  the  presence  of  V1,  a,  and  0  in  the  free  energy  function  allows  for 
contribution  of  surface  or  interfacial  energy  of  grain  boundaries,  for  example,  in 
manifesting  certain  length  scale  (surface  area  to  bulk  volume)  effects  on  flow  stress, 
provided  the  dislocation-disclination  content  of  such  boundaries  is  described.  Such 
effects  are  dominant  for  nanocrystalline  materials. 

The  dissipation  inequality  (44)  may  be  re-written  as  follows,  upon  introduction  of 
generalized  second-rank  thermodynamic  forces  a  and  Jv. 


<A  <  LapMT!  -  *>>~U  ~  ^  (9“r,a)  -  NT  +  ojZ  +  t'X’ 


(46) 


where  the  final  two  terms  vanish  upon  consideration  of  the  material  time  derivative 
of  Eqs.  (38).  Substituting  (45)  into  (46)  and  assuming  that  (46)  must  hold  as  C  and  T 
are  varied  independently  of  the  other  constitutive  variables  (e.g.,  a  reversible  thermo¬ 
elastic  process),  we  arrive  at  the  usual  constitutive  relations,  i.e., 


0C 


=  .J'F 


<e—  1  a  ^abj^e—  1  /? 


N=- 


a/? 


01 } 
ST’ 


(47) 


with  S  a  macroscopic  elastic  second  Piola-Kirchhoff  stress,  as  well  as  the  strong  form 
of  the  reduced  dissipation  inequality, 


(p,Lp)  +  (pv^y)  +  (vpv^w1)  +  +  (fX) 


(48) 


where  the  spin  w'  and  plastic  velocity  gradient  LP  referred  to  configuration  b  are 
found  from 
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f'F'1  =  VV'~l  +  VrJr^V1-1,  Lp  =  F'FPF|,  lFi  (49) 

w1 

and  where  the  Eshelby-type  stress  (cf.  Eshelby,  1951,  1975;  Maugin,  1994;  Le  and 
Stumpf,  1996c), 

P  =  MT  -  =  Jl Te-'  (Eg  -  t/H)  Fe  =  J  'f  (FtT  G  -  iA0l0)  F“‘  (50) 


is  defined  as  the  push-forward  of  the  mixed-variant  reference  configuration  energy- 
momentum  tensor  P,  or  equivalently  as  the  pull-back  of  the  current  configuration 
Eshelby  stress  tensor  p.  In  Eq.  (50),  10,  1,  and  1  denote  mixed-variant  identity  maps 
on  configurations  b0,  b,  and  b,  respectively,  T  =  .S'1  is  the  usual  first  Piola-Kirchhoff 

~  T 

stress,  M  is  the  transposed  Mandel  stress,  g  is  the  current  configuration  metric  ten¬ 
sor,  G  is  the  reference  configuration  metric  tensor,  and  i/r  is  the  free  energy  measured 
per  unit  current  volume.  We  enforce  the  notion  that  the  geometrically  necessary  de¬ 
fect  density  tensors  a  and  6  contribute  neither  storage  or  dissipation  on  the  right  side 
of  (48),  resulting  in 


_  0i }  _  _  0tA 

“as’  M _00 


and  leading  to 


(51) 


Dissipation  from 
plastic  defect  fluxes 
(e.g.  dislocation  glide) 


Dissipation  from 
microscopic  heterogeneity 
(e.g.  dislocation  clustering) 


/viPVi-\Wi\  + 

S - V - " 

Dissipation  from 
microscopic  heterogeneity 
(e.g.  disclination  spin) 


Dissipation  from 

geometrically  necessary  defect  densities 


Energy  storage  from 
microscopic  heterogeneity 


+ 


Energy  storage  from 
statistically  stored  defect  densities 


Choice  of 

intermediate  metric  tensor 


(52) 


We  have  described  the  physics  behind  each  term  in  (52)  apart  from  the  final  term  re¬ 
lated  to  g.  If  we  make  the  (usual)  constitutive  assumption  ga/3  =  <5a/J,  where  is  the 
covariant  identity  map,  this  term  vanishes.  As  discussed  by  Clayton  et  al.  (2004b), 
other  choices  for  gx/j  may  be  advantageous  for  describing  situations  in  which  internal 
stresses  associated  with  defect  arrangements  (e.g.,  at  flexing  subgrain  walls)  may 
arise  from  or  be  amplified  by  applied  loads  (Gibeling  and  Nix,  1980;  Argon  and 
Takeuchi,  1981;  Kassner  et  al.,  2002).  Additionally,  note  that  from  (38)  we  have 
(ffT,ap)  =  —  (<fT,a)  and  (/iT,0p)  =  — (/iT,0),  meaning  that  the  final  two  terms  on 
the  left-hand  side  of  inequality  (52)  may  be  replaced  by  terms  denoting  energy  stor¬ 
age  associated  with  geometrically  necessary  defects  (see  later  Eq.  (90)),  as  a  matter  of 
preference  in  addressing  the  issue  of  energy  storage  versus  dissipation  as  a  function 
of  scale  (cf.  Wei  and  Hutchinson,  1999). 

The  general  philosophy  adopted  here  is  prescription  of  balance  laws  for  the 
energetic  forces  I  (macroscopic  Cauchy  stress),  a  (microscopic  force  conjugate  to 
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geometrically  necessary  dislocations  within  the  volume  element),  and  Ji  (microscopic 
force  conjugate  to  geometrically  necessary  disclinations  within  the  volume  element) 
as  well  as  the  temperature  T.  Additional  evolution  equations  are  then  required  to 
specify  the  time  history  of  dissipative  kinematic  variables  F1  and  Fp  and  scalar  inter¬ 
nal  variables  ep  and  e,.  Force  conjugates  J^j,  Sf,  and  ^  associated  in  (52)  with  energy 
storage  (i.e.,  residual  energy  from  cold  working)  are  not  required  to  satisfy  any  addi¬ 
tional  balance  laws. 

From  (39)  to  (41),  the  macroscopic  Cauchy  stress  I  is  supposed  to  obey  the  stan¬ 
dard  linear  and  angular  momentum  balances  and  reflects  the  average  traction  carried 
by  a  local  crystalline  volume  element  in  the  current  configuration.  The  microforces  a 
and  ji  reflect  higher  order  moments  of  the  microscopic  traction  distribution  sup¬ 
ported  by  the  volume  element  (see  e.g.,  Kroner,  1963)  and  do  not  explicitly  enter 
the  macroscopic  momentum  (41)  or  energy  (42)  balances  in  our  theory.  Instead,  con- 
travariant  mixed-configurational  versions  of  these  forces, 

=  Je-lF-arx^sgp5F%  (53) 

satisfy  coupled  microscopic  momentum  balances  independent  from  (40).  The  higher 
order  microscopic  balance  relations  are  intended  to  apply  to  subvolumes  within  the 
considered  macroelement  pertinent  to  the  macroscopic  momentum  balance  relations 
in  (41).  At  these  finer  scales,  heterogeneous  features  such  as  grains,  second  phases, 
subgrains,  or  other  dislocation  substructures  are  expected  to  engender  gradients  of 
a  and  6.  This  clarification  of  macro-  and  micro-effectively  parallels  the  kinematical 
argument  of  Section  2  in  that  F'  reflects  heterogeneity  at  subscales  within  the  consid¬ 
ered  volume  element  subjected  to  the  multiplicative  deformation  gradient  decompo¬ 
sition  F  =  FeF‘Fp  =  F^FP. 

Please  note  that  a  is  a  first-order  moment  stress  (i.e.,  couple  stress)  with  units  of 
Force! Length,  whereas  /i  is  a  second-order  moment  stress  (i.e.,  hyperstress)  with  units 
of  Force.  Extending  standard  stress  tetrahedron  arguments  for  micropolar  media 
with  higher-order  stresses  (cf.  Malvern,  1969),  we  postulate  the  following  global 
angular  momentum  balances,  restricted  here  to  the  microscopically  quasi-static  case 
(i.e.,  no  micro-inertia)  and  in  the  absence  of  microscopic  body  forces: 


/  ta  da  =  0,  /  (s  :  (x  <g>  tff)  +  t,,)  d a  =  0,  (54) 

J  a  J  a 

where  t„  =  (a,  n)  and  =  (ju,  n)  are,  respectively,  first-  and  second-order  couple  trac¬ 
tion  vectors  with  components  referred  to  b  and  acting  on  spatial  surface  a  with  out¬ 
ward  normal  n,  and  x  =  Vv'  'x  are  spatial  coordinates  pulled  back  to  b  under  the 
Cauchy-Born  approximation.  Localized  forms  of  (54)  then  follow  naturally  from 
application  of  the  divergence  theorem  in  the  current  configuration: 

<  =  0,  e^Ff-1Vi  +  4  =  0,  (55) 

where  /C  =  and  vertical  bars  denote  covariant  differentiation  with  respect  to 

spatial  coordinates  x.  Relations  (55)  are  analogous  to  the  micropolar  elastic  balance 
laws  used  in  the  finite  strain  theory  of  Minagawa  (1979)  and  the  linearized  theory  of 
Eringen  and  Claus  (1970).  We  note  that  if  \Jt  is  quadratic  in  each  of  os  and  0,  then  a 
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and  /i  in  (55)  are,  respectively,  linear  in  configurational  mappings  of  a  and  0,  and  the 
so-called  “microforce  balances”  in  (55)  are  then  interpreted  in  terms  of  covariant 
spatial  gradients  of  geometrically  necessary  dislocations  (55)i,  pertaining  to  pile- 
ups  and  dipole  structures,  and  disclinations  (55)2,  contributing  to  nonsymmetric  mo¬ 
ment  stress  <r. 

The  usual  heat  equation  is  suggested  for  the  temperature  field: 

qA  =  -kT\BGl!A,  (56) 


where  we  have  assumed  isotropic  conduction  (a  typical  assumption  easily  extended 
to  the  anisotropic  case),  with  thermal  conductivity  k,  posed  above  in  the  reference 
coordinate  frame. 

The  time  rate  of  Fp  is  described  by  Eq.  (9),  i.e.,  FP  =  LPFP,  which  in  turn  requires 
constitutive  equations  for  slip  rates  yl  of  (10).  In  rate  dependent  crystals,  slip  is  typ¬ 
ically  dictated  by  a  viscoplastic  flow  rule,  in  our  case  of  the  following  form  in  con¬ 
figuration  b  : 


\t),  (57) 

where  the  driving  force  for  slip,  p'  =  rPsF'~la-P  pF1^,  is  a  resolved  scalar  Eshelby 
stress  on  system  i  such  that 


P,F' 


m 


F1-'  )  = 


(58) 


is  the  contribution  of  the  plastic  shearing  rates  to  the  dissipation  inequality  (52).  The 
variables  p‘s  and  Tfs  are  statistically  stored  defect  densities  assigned  to  the  1th  slip  sys¬ 
tem,  i.e.,  line  lengths  per  unit  volume  in  configuration  b,  such  that 

J*Ps  =  Ps  =  J2  Ps’  Jif!s  =  Vs  =  Y  (59) 


with  J'  the  Jacobian  determinant  of  F1.  In  (57),  statistically  stored  defects  generally 
contribute  to  self-hardening  ( i  =  j)  and  latent-hardening  (i  p  j).  Geometrically  neces¬ 
sary  defect  densities  also  contribute  to  hardening  on  each  system,  and  as  a  result, 
their  counterparts  in  configuration  b,  defined  as 


9SP  =  J'F‘: 


V'-1? 


(60) 


are  included  in  the  flow  rule  (57).  Please  note  that  the  net  defect  density  tensors  in 
(60)  are  not  to  be  confused  with  the  total  defect  density  tensors  at1,  and  0I  introduced 
in  Eq.  (8).  It  is  also  noted  that  within  our  framework,  the  plastic  velocity  gradient 
may  alternatively  be  described  by  evolution  laws  for  defect  fluxes  and  generation 
rates  of  Eqs.  (7)  or  (10)  (cf.  Acharya  (2001,  2003))  rather  than  evolution  equations 
for  slip  rates  as  in  (57).  The  former  approach  appears  more  conducive  to  correlation 
with  numerical  experiments  based  on  discrete  dislocation  dynamics  (cf.  Zbib  and 
De  La  Rubia,  2002)  or  atomic-scale  interactions  (cf.  Nazarov  et  al.,  2000),  while 
the  latter  approach  appears  more  conducive  to  correlation  with  macroscopic 
experiments. 
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Specific  guidelines  describing  how  dislocation  densities  may  participate  in  the  flow 
rule  (i.e.,  strain  hardening)  can  be  found  in  previous  works  by  the  co-authors 
(Bammann,  2001;  Regueiro  et  al.,  2002;  Clayton  et  al.,  2004a).  Disclination  densities 
(specifically,  partial  disclination  dipoles)  may  contribute  to  hardening  and/or  geomet¬ 
rical  softening  depending  upon  the  particular  material  under  consideration  and  its 
deformation  history,  as  discussed  by  Pccherski  (1983, 1985)  and  Seefeldt  and  co-work¬ 
ers  (Seefeldt  and  Klimanek,  1997,  1998;  Seefeldt,  2001).  The  right  stretch  tensor  U1  of 
F1  (i.e.,  the  counterpart  of  left  stretch  V'  in  configuration  b)  is  included  in  (57)  to 
account  for  contributions  of  heterogeneity  of  inelastic  deformation  within  the  local 
volume  not  engendered  by  that  volume’s  (average)  statistically  stored  and  geometri¬ 
cally  necessary  defect  densities,  such  as  higher  order  moments  of  dislocation  densities 
(Kroner,  1973;  Hartley,  1975)  contained  within  the  element.  For  example,  a  backstress 
may  arise  in  conjunction  with  piled-up  dislocations  at  a  misoriented  subgrain  bound¬ 
ary.  The  corresponding  inhibiting  effect  on  plastic  flow  would  then  manifest  in  (57) 
through  the  inclusion  of  6  and  U1,  the  former  reflecting  the  misorientation  boundary 
itself  and  the  latter  accounting  for  a  spatial  gradient  of  dislocation  density  in  the  pile- 
up  direction.  Our  treatment  could  be  further  extended  to  explicitly  include  higher  or¬ 
der  gradients  or  moments  of  defect  densities  in  the  flow  rule  or  yield  conditions,  albeit 
at  the  expense  of  increased  model  complexity.  Also,  when  inelastic  volume  changes  are 
of  interest  (e.g.,  high  temperature  creep  conditions  characterized  by  vacancy  migra¬ 
tion),  a  nonvanishing  non-slip  term  should  be  included  in  kinematic  relation  (9), 
with  its  evolution  exhibiting  the  same  general  constitutive  dependency  as  f  in  (57). 

Since  statistically  stored  defects  accumulate  in  response  to  the  motion  of  defects 
and  their  interactions,  their  evolution  equations  are  assigned  the  same  general  con¬ 
stitutive  dependencies  as  the  shearing  rates,  i.e., 


(61) 


Kinetic  relations  such  as  (61)  are  strongly  material  dependent.  More  specific  equa¬ 
tions  dealing  with  coupled  evolution  of  dislocation  and  disclination  densities  -  which 
may  be  considered  as  particular  implementations  of  the  general  relations  (61)  -  may 
be  found  in  the  cited  references  (Pccherski,  1983,  1985;  Seefeldt  and  Klimanek,  1997, 
1998;  Seefeldt,  2001).  An  analogous  general  prescription  is  given  for  the  time  rate  of 
F1,  referred  to  configuration  h\ 


(62) 


where  P^  =  F~UP  pF'?  is  the  Eshelby  stress  tensor  on  b.  As  introduced  in  Eq.  (12),  F1 
accounts  for  lattice  heterogeneity  within  the  crystalline  volume  element  and  also  dic¬ 
tates  the  rate  of  rotation  of  the  slip  planes  and  directions  arising  from  the  generation 
and  flux  of  mobile  disclinations  traversing  the  element,  and  as  a  result,  will  account 
for  the  diffusion  of  crystallographic  texture  observed  as  subgrain  cells  -  delineated  in 
our  present  framework  by  walls  of  partial  disclination  dipoles  -  emerge  and  then  ro¬ 
tate  relative  to  one  another  (see  e.g.,  Butler  and  McDowell,  1998;  Peeters  et  al.,  2001; 
Hughes  et  al.,  2003).  From  Fig.  4  and  the  definition  F :/’  =  FeF‘,  the  rotational  part  Re 
of  Fe  depicts  the  average  rotation  undergone  by  the  lattice  between  unloaded  state  b 
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and  the  current  configuration  b,  while  the  rotational  part  R1  of  F1  accounts  for  addi¬ 
tional  rotation  of  the  lattice  due  to  defect  fields  embedded  within  the  material  in  re¬ 
laxed  configuration  b.  Notice  also  that  F1  affects  the  slip  system  orientations  through 
its  inclusion  in  the  total  lattice  deformation  Fy\  since  the  slip  directions  s'  and  slip 
plane  normals  m'  are  pushed  forward  to  current  configuration  b  from  intermediate 
configuration  b  as 

s'  =  F^s',  m!  =  m'F^1.  (63) 

The  specific  form  of  the  constitutive  equation  for  £'  for  a  given  material  and  scale  of 
resolution  (e.g.,  /rcf  in  (1))  depends  upon  the  particular  arrangement  of  defects  with¬ 
in  the  crystal  and  how  their  local  lattice  strain  fields  interact  and  contribute  to  the 
motion  of  the  element’s  external  boundary.  Discrete  dislocation  simulations  or  atom¬ 
istic  studies  conducted  in  a  manner  analogous  to  the  continuum  simulations  de¬ 
scribed  in  Clayton  and  McDowell  (2003)  are  expected  to  provide  more  insight 
into  developing  Eq.  (62)  for  particular  materials. 

The  micropolar  kinematic  variable  Q  describes  spatial  gradients  of  lattice  rotation 
within  the  crystalline  volume  element  arising  from  distributed  disclinations.  Consider 
parallel  transport  of  the  lattice  director  d„,  over  the  small  distance  dx  to  a  new  ori¬ 
entation  <f ,  conducted  with  respect  to  the  covariant  derivative  (19).  To  isolate  the 
effects  of  Q,  we  restrict  attention  now  to  the  case  when  Fy  «  1  and  spatial  gradients 
of  the  lattice  deformation  Fy  are  small.  Applying  the  summation  convention  over 
covariant  Cartesian  coordinates,  we  then  may  write 

da  da  +  &a,bdxb  d„  -\-  Qbacdxb dL.  ^  for  +  Qh\ar\dxj)  d,  (pacdc .  ( 64 ) 

where  Qhac  -  Qb[ac]  from  Eq.  (21),  and  the  rotation  matrix  <Pac  satisfies  <FT  ~  <I>  1 
and  det(<D) «  1  for  small  magnitudes  of  Qbacdxh,  i.e.,  small  relative  lattice  rotations 
such  as  those  occurring  across  low-angle  grain  boundaries.  We  see  from  (64)  that 
Qbac  acts  as  an  effective  “gradient”  of  rotation  in  the  spatial  direction  xh.  However, 
Q  is  not  a  true  spatial  gradient  since  Qb\ac\d.-Q<\ac\b  fo  0  ->  Qh[ac]  fo  i)[aclb,  where  i)[ac] 
is  a  skew  matrix  that  can  exist  only  when  the  compatibility  conditions  Qb[ac],d- 
Qd[ac],b  are  met.  From  Eq.  (26)  and  the  second  of  (29),  enforcement  of  such  condi¬ 
tions  would  in  fact  yield  a  null  disclination  density  0  in  the  context  of  linear  approx¬ 
imation  (64).  Note  that  the  nine  independent  degrees-of-freedom  of  Q  could  be 
determined,  in  principle,  from  knowledge  of  incremental  changes  in  lattice  rotation 
due  to  disclinations,  based  upon  measured  gradients  of  intragranular  misorientation, 
sometimes  referred  to  as  “disorientation”  (Barton  and  Dawson,  2001;  Bergugnat, 
2002;  Kassner  et  al.,  2002). 

Our  framework  is  complete  upon  introduction  of  explicit  relationships  between 
the  macroscopic  rotation  attributed  to  defects,  R1,  and  the  microscopic  defect  density 
measures,  e.g., 

Ri  =  Ri(a,0,8p,e,).  (65) 

Please  notice  that  (65)  is  introduced  independently  of  kinetic  relation  (62),  the  latter 
implicitly  containing  an  evolution  equation  for  the  spin  due  to  the  time  rate  of  R', 
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i.e.,  W'  of  (49),  which  in  turn  contributes  to  the  dissipation  rate  in  (52)  but  not  to 
energy  storage  in  i }  of  (45).  Also  note  that  no  corresponding  kinematic  equation 
analogous  to  (65)  is  given  for  the  stretch  V1,  which  is  included  explicitly  in  the  free 
energy  function  of  (45)  to  account  for  higher  order  moments  of  distirubutions  of 
lattice  defects  that  lead  to  internal  heterogeneity  of  lattice  deformation  within 
the  local  volume  element,  such  as  dislocation  pile-ups,  for  example,  that  are  not 
accounted  for  by  variables  a,  0.  ep,  and  en.  A  variant  of  (65)  follows  from  Lardner 
(1973,  1974),  if  we  assume  that  R1  is  determined  by  the  change  in  orientation  of  a 
lattice  director  vector  upon  parallel  transport  around  a  closed  circuit  encircling 
area  a: 

R‘$  =  SI  (dab  +  J  ^  R‘bcdgCdene  da)  <5*  =  d«  +  5 ^  f  C^~ubebe^hnh  da 

« a?  +  (66) 

where  S]  and  dhp  are  shifters  between  coincident  Cartesian  spatial  and  intermediate 
frames,  and  where  in  the  final  equality  of  (66),  we  have  assumed  a  constant  lattice 
stretch  over  a,  such  that  the  net  Frank  vector  Qs  =  \a68hnh  da  plays  the  role  of  the 
axial  vector  corresponding  to  the  rotation  tensor  R1  in  (66),  rigorously  orthogonal 
only  in  the  limit  of  small  rotations.  Notice  that  definition  (66)  is  non-unique  in 
the  sense  that  R‘  depends  upon  the  choice  of  oriented  area  element  nda.  In  our  crystal 
plasticity  framework,  it  is  natural  to  generalize  the  second  of  (66)  to  a  discrete  aver¬ 
age  taken  over  potentially  active  slip  planes,  as  R'  contributes  to  the  average  rotation 
of  the  lattice  directors  and  the  slip  planes  (e.g.,  Eqs.  (15)  and  (63)).  Such  an  approach 
gives 

0*V)J,  (6V) 

where  k  is  the  number  of  slip  systems  and  the  factor  f2re[  corresponds  to  the  effective 
area  over  which  the  disclinations  act.  In  the  limit  that  our  local  crystal  volume  ele¬ 
ment  shrinks  to  infinitesimal  size,  Aef  — >  0,  and  effects  of  disclinations  can  no  longer 
be  resolved  by  the  deformation  gradient  field.  Particular  Eq.  (67)  is  deemed  an 
appropriate  approximation  of  (65)  under  the  following  conditions:  (i)  the  elastic  lat¬ 
tice  stretch  and  net  Frank  vector  are  small,  (ii)  disclinations  dominate  the  net  resid¬ 
ual  lattice  rotation  relative  to  geometrically  necessary  dislocations,  and  (iii)  the  net 
influences  on  lattice  rotation  from  statistically  stored  defects  contained  within  the 
volume  element  cancel  out.  Furthermore,  when  the  residual  lattice  stretch  associated 
with  microstructural  heterogeneity  (V1)  is  negligible,  (12)  reduces  to  the  kinematic 
description  of  Lardner  (1973,  1974)  and  Pgcherski  (1983,  1985),  i.e.,  F  =  FeR‘Fp. 
On  the  other  hand,  in  the  most  general  case,  the  multiplicative  decomposition  (12) 
may  be  extended  to  F  =  (Fe)(F'(1)F'(2) . . . F^^F^F^  . . .  Fp^)  where  the  residual  lat¬ 
tice  deformation  F‘  and  lattice-preserving  plastic  deformation  Fp  are  each  resolved 
into  components  associated  with  k  distinct  classes  of  defects  (e.g.,  statistically  stored 
dislocations,  geometrically  necessary  dislocations,  disclinations,  etc.).  Then  from  (45) 
and  (52),  each  F'(t)  would  contribute  to  free  energy  storage  (through  a  conjugate 


R*  =  51  +  5y 
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thermodynamic  force)  and  dissipation  (through  the  macroscopic  stress  power),  while 
each  would  only  contribute  to  dissipation. 

4.  Model  summary  and  elements  of  implementation 

Summarized  below  are  the  primary  relationships  comprising  the  theory,  assuming 
for  simplicity  isochoric  plastic  deformation  Lp  characterized  by  slip-type  kinematics 


and  disclination-dominated  residual  lattice  rotation,  R1: 

Macroscopic  kinematics 

(deformation  gradient)  F  =  FeF'Fp  =  F^FP,  (68) 

(plastic  velocity  gradient)  Lp  =  FPFpl  =  y^y's'  <g>  m',  (69) 

i 

(slip  system  directors)  s' =  F^s',  m'  =m'¥se~l,  (70) 

(total  lattice  strain)  Cfb  =  (71) 

Microscopic  kinematics 

(linear  connection)  Pfh  =  F*aF%~Vi  +  Qcab,  (72) 

(micro  -  rotation)  Qcba  =  Qc[ba]  =  QjC%,  (73) 

(torsion/dislocations)  T;ah  =  Tgb]  =  Pcb  -  Pbc  =  scbdotad,  (74) 

(curvature/ disclinations)  Rabcd  =  R[ab][cdl  =  2V[cQd][ba]  +  TfdQe[ba] 

=  ZbaePdfOef-  (75) 

Micro-macroscopic  kinematics 


/A  _  , 

(disclination  lattice  rotation)  R =  <5|  +  baadepCse~lab sbeg  ( 


Macroscopic  balance  laws 

(76) 

(mass)  p0  =  pJ, 

(77) 

(linear  momentum)  Iah  ,b+  pfa  =  pxa. 

(78) 

(angular  momentum)  sabcShc  =  0, 

(79) 

(traction  boundary  conditions)  f  =  Iabnb. 

(80) 
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Microscopic  balance  laws 

(first  order  micro  —  momentum)  a 


.5*  =  Q 


(second  order  micro  -  momentum)  ^a/b  +  ptb  =  0, 

(traction  boundary  conditions)  (7ff)a  =  <raA«i,  (tbY  =  pahnb. 


Thermodynamics 

(free  energy)  i }  =  \j/(c*,  Wl,i,e,sp,s^T,gj, 

(internal  energy  rate)  U0  =  SAagabFhA  -qfA+H0, 

(heat  flux)  qA  =  -kT\BGBA, 

(macrostress)  ~S*P  =  2  =  fF^YahFe-'s\ 

dib 

(entropy)  N  =  -  —  , 

/  •  .  _  01 }  _  0i } 

(microstresses)  <r  =  — ,  p  = 

(dissipation  inequality)  ^P,LP^  +  ^PV'  V*^  +  ^V'PV'  ^W* 


(81) 

(82) 

(83) 

(84) 

(85) 

(86) 

(87) 

(88) 
(89) 


(90) 


Kinetics 

(plastic  slip  rates)  f  =  y'(pl,  p's,  rfs ,  a,  6 ,  U\  T) ,  (91) 

(micro  -  heterogeneity)  F'~‘f'  =  £'  =  L'(P,  a,  0,  U‘,  T),  (92) 

(statistically  stored  dislocations)  p's  =  ps(p\  p's ,  ff's,  a,  0,  U1,  T) ,  (93) 

(statistically  stored  disclinations)  rj's  =  rj'siP'i  Psi  Ys-  *  ^  U1,  71) .  (94) 


Some  general  comments  are  now  given  regarding  aspects  of  our  constitutive  the¬ 
ory  that  differ  from  others  in  the  literature.  Our  macroscopic  ((78)  (80))  and  micro¬ 
scopic  ((8 1)  (83))  force  balance  equations  are  independent,  with  the  symmetry  of  the 
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stress  tensor  (79)  and  usual  structure  for  internal  energy  rate  (85)  maintained.  This  is 
in  contrast  to  the  couple  stress  theory  conceived  by  Fleck,  Hutchinson,  and  co-work¬ 
ers  (Fleck  and  Hutchinson,  1993;  Fleck  et  al.,  1994;  Gao  et  al.,  1999)  that  permits  a 
nonsymmetric  Cauchy  stress.  Gurtin  (2002,  2004)  developed  a  gradient  crystal  plas¬ 
ticity  theory  preserving  the  symmetry  conditions  on  the  macroscopic  stress  and 
simultaneously  featuring  microforce  balances  on  discrete  slip  systems  whereby  geo¬ 
metrically  necessary  dislocations  (in  terms  of  both  first  and  second-order  spatial  gra¬ 
dients  of  a  lattice  deformation  map  Fe  1 )  may  naturally  contribute  to  a  backstress  on 
each  system.  In  Gurtin’s  theory,  slip  rates  perform  mechanical  work  on  the  boundary 
of  the  body  and  explicitly  enter  the  global  energy  balance  (virtual  power  principle) 
from  which  governing  equations  are  derived.  In  our  theory,  the  micro-variables  do 
no  mechanical  work  (they  are  absent  in  (85)). 

In  the  present  work  we  have  postulated  the  existence  of  microscopic  force  balances 
from  the  outset,  following  consideration  of  balances  of  higher  order  moments  (54). 
Alternative  approaches  for  constructing  balance  relations  in  higher  order  theories 
often  employ  variational  or  virtual  work  principles  (e.g.,  Toupin,  1964;  Teodosiu, 
1967;  Le  and  Stumpf,  1996a, c;  Gurtin,  2002).  Our  thermodynamic  force  relations 
(89)  emerge  upon  division  of  the  total  defect  density  tensors  into  elastic  (storage) 
and  plastic  (dissipative)  parts,  as  in  Eq.  (38).  Le  and  Stumpf  (1996a,c)  made  an  alter¬ 
native  prescription,  dividing  the  conjugate  defect  forces  into  dissipative  and  energetic 
parts,  the  former  associated  with  evolution  laws  for  defect  fluxes  and  the  latter  with 
stress-type  balance  laws  for  recoverable  elastic  deformations  and  defect  densities. 

Prior  to  implementation  and  subsequent  solution  of  initial-boundary  value  prob¬ 
lems,  specific  forms  of  the  free  energy  function  (84)  and  kinetic  relations  (91)  (94) 
must  be  selected  for  the  material  under  consideration.  The  framework  presented  here 
establishes  the  need,  and  serves  as  a  guide,  for  developing  appropriate  experiments, 
both  physical  and  numerical,  necessary  for  simulation  of  a  generalized  continuum 
accounting  for  the  generation,  motion,  and  interaction  of  dislocation  and  disclina- 
tion  defects.  The  dependence  of  the  free  energy  upon  periodic  defect  distributions 
can  be  extracted  from  homogenization  of  atomistic  solutions,  for  example  (Chung, 
2004).  A  non-convex  potential  (84)  will  promote  the  development  of  heterogeneous 
deformation  fields  in  a  body  subjected  to  nominally  homogeneous  initial  and  bound¬ 
ary  conditions,  as  discussed  by  example  in  Appendix  A.  Note  that  normalizing 
length  scale  parameters  or  constants  (Gao  et  al.,  1999;  Regueiro  et  al.,  2002;  Al- 
Rub  and  Voyiadjis,  2004;  Huang  et  al.,  2004)  associated  with  the  dislocation  and  dis- 
clination  density  tensors  are  expected  to  enter  the  free  energy  and  kinetic  relations,  as 
these  defect  densities  are  not  dimensionless. 

From  the  perspective  of  numerical  solution  of  boundary  value  problems,  Eqs. 
(68)-(94)  are  formidable.  Consider  this  framework  in  comparison  to  that  of  standard 
crystal  plasticity  (e.g.,  Asaro,  1983).  The  additional  kinematic  unknowns  present  in 
the  former  relative  to  the  latter  are  F1  (  Eq.  (68))  and  Q  (Eq.  (73)),  each  having  nine 
independent  components.  Note  that  the  defect  density  tensors  can  be  derived  directly 
from  Fy  and  Q  via  (74)  and  (75).  The  corresponding  additional  equations  in  our 
framework  are  kinematic  relations  (76)  (3  equations  relating  the  residual  lattice 
rotation  and  the  disclination  density,  which  in  turn  is  a  function  of  Q  and  its  spatial 
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gradients),  microforce  balances  (81)  and  (82)  (a  total  of  6  equations),  and  the  kinetic 
equations  for  F'  in  (92)  (9  equations).  When  geometrically  necessary  disclinations  are 
absent,  Q  =  0,  R1  =  1,  geometrically  necessary  dislocations  are  found  from  curl  of  the 
(inverse)  lattice  deformation  gradient  via  (72)-(74),  and  satisfaction  of  the  micro¬ 
force  balances  is  trivially  ensured  upon  proper  choice  of  free  energy  dependency 
on  dislocations  (e.g.,  a  dislocation  driving  force  a  proportional  to  a  curl  of  a  tensor 
field  will  be  divergence-free).  For  example,  a  particularly  simple,  yet  illustrative, 
form  for  the  free  energy  function  (84)  in  the  absence  of  disclination  densities  is  the 
following: 

-A  =  {c:p-~g^c^(^xS-~gxS)  + 

+  ^2  (zp)  +  ^3(^)1  (95) 

where  Cafo<'  are  temperature-dependent  elastic  moduli,  and  where  t/q,  1 ]/2,  and  t/q  are 
functions  of  their  arguments,  with  1 /q  accounting  for  the  specific  heat  capacity  of  the 
material.  Notice  that  the  second  term  in  (95)  consists  of  a  quadratic  form  of  the  sym¬ 
metric  part  of  the  dislocation  density  tensor  5,  pre-multiplied  by  an  effective  elastic 
shear  modulus  Ji  and  normalized  by  squared  length  parameter  {.  In  this  case,  setting 
the  micro-traction  (hyperstress)  boundary  condition  =  0  for  the  second  of  (83)  in 
conjunction  with  the  prescription  IF  1 R  =  0  for  spin  components  of  rate  Eq.  (92) 
leads  to  the  (correct)  trivial  solution  Q  =  0  in  (72)  and  R1  =  1  in  (76),  assuming  for 
initial  conditions  a  defect-free  reference  lattice.  In  other  words,  substituting  (95)  into 
(81)  and  (82)  yields  the  identically  satisfied  expressions  n®  =  l^2£bcdF^~^  =  0  and 
EapIF^-ip<jxb  =  Je~1p/2eafiF,~lia,^F‘^  =  0,  respectively,  written  here  in  Cartesian 
spatial  coordinates.  Notice  that  the  defect  densities  (i.e.,  higher  order  deformation 
gradients)  can  enter  the  kinetic  relations  (91)— (94),  effectively  rendering  these  partial 
differential  equations  in  time  and  space.  Further  elaboration  upon  the  role  of  higher- 
order  traction  boundary  conditions  (83)  and  their  kinematic  counterparts  is  deferred 
to  future  work  wherein  specific  free  energy  and  kinetic  relations  will  be  developed  for 
plastically  deforming  crystals  containing  evolving  distributions  of  dislocations  and 
disclinations. 


5.  Concluding  remarks 

A  general  finite  strain  theory  incorporating  dislocation  and  disclination  defects 
has  been  introduced,  based  primarily  on  volume  averaging  considerations.  Novel 
features  of  the  kinematic  description  include  a  three-term  decomposition  of  the  aver¬ 
age  deformation  gradient  for  a  crystalline  volume  element  and  an  independent 
micropolar  rotation  variable  associated  with  incompatible  spatial  gradients  of  lattice 
rotation  within  the  element.  The  latter  may  be  attributed  to  misoriented  subgrain 
cells  separated  by  walls  of  partial  disclination  dipoles.  The  micropolar  rotations 
execute  independently  of  the  continuum  rotations  associated  with  elastic  and  inelas¬ 
tic  components  of  the  deformation  gradient,  and  are  intended  to  capture  physics 
that  geometrically  necessary  dislocations  associated  with  first  order  gradients  of 
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incompatible  elastic  or  plastic  deformations  cannot  represent  alone.  For  example, 
these  micro-rotations  may  evolve  even  under  nominally  homogeneous  applied  defor¬ 
mations  to  single  crystals,  in  contrast  to  geometrically  necessary  dislocation  theories 
requiring  heterogeneity  of  deformation  at  the  scale  of  the  crystal  to  drive  the  evolution 
of  defects  (see  the  discussion  in  Ortiz  and  Repetto,  1999;  Ortiz  et  ah,  2000,  whose  lam¬ 
ination  model  also  intends  to  represent  such  effects).  Our  theory  separates  the  role  of 
incompatibility  into  contributions  from  higher  order  deformation  gradients  (incor¬ 
porated  primarily  in  the  torsion  or  geometrically  necessary  dislocation  tensor)  and 
subgranular  lattice  rotation  gradients  (incorporated  in  the  curvature  or  geometri¬ 
cally  necessary  disclination  tensor),  although  the  constitutive  equations  ultimately 
dictating  their  evolution  are  in  general  strongly  coupled.  The  kinematic  framework 
has  strong  ties  to  geometrically  oriented  defect  field  theories  of  the  mid-  to  late- 
twentieth  century  (Bilby  et  al.,  1955;  Noll,  1967;  Lardner,  1973;  Minagawa,  1979; 
De  Wit,  1981),  and  reduces  to  a  formulation  consistent  with  previous  gradient-based 
dislocation  theories  of  the  co-authors  (cf.  Bammann,  2001;  Regueiro  et  al.,  2002; 
Clayton  et  ah,  2004a)  when  micropolar  rotations  and  disclinations  are  excluded. 
Our  theory  admits  the  coupling  of  applied  and  residual  elastic  strain  energy  densities, 
arising  physically  from  the  bowing  of  dislocation  structures  due  to  applied  loads, 
through  the  appropriate  choice  of  free  energy  function  (45).  Higher  order  micro¬ 
stresses  satisfy  balance  laws  independent  of  those  associated  with  the  macroscopic 
Cauchy  stress,  and  exhibit  a  form  comparable  to  classical  couple  stress  theory 
(Eringen  and  Claus,  1970;  Minagawa,  1979).  Kinetic  relations  are  left  in  a  general 
format  applicable  to  generic  metallic  crystals.  We  expect  the  additional  assertion 
of  maximal  dissipation  rate  to  serve  as  an  heuristic  tool  for  partitioning  disloca¬ 
tion  and  disclination  kinetics.  Also  required  is  a  more  thorough  examination  of 
the  supplemental  boundary  conditions  (83)  associated  with  the  lattice  arrangement 
(cf.  Acharya  et  ah,  2004). 

The  model  is  expected  to  capture  grain  size  effects  on  initial  yielding  and  strain 
hardening,  the  former  via  proper  choices  of  initial  conditions  for  geometrically  nec¬ 
essary  dislocations  and  disclinations  and  their  inclusion  in  the  initial  slip  resistances, 
the  latter  via  evolution  of  the  crystal  defects  (cf.  Evers  et  ah,  2002)  and  their  incor¬ 
poration  in  the  viscoplastic  flow  rule.  Since  relatively  strong  lattice  rotation  gradients 
are  expected  in  the  vicinity  of  grain  boundary  regions  (Bergugnat,  2002)  and  may  be 
correlated  to  measured  hardness  increases  near  such  boundaries  (Soifer  et  ah,  2002), 
initial  conditions  on  the  micro-rotation  variable  and  defect  densities  should  reflect 
larger  initial  yield  stresses  in  relatively  fine-grained  microstructures  that  exhibit  more 
grain  boundary  surface  area  per  unit  volume  than  relatively  coarse-grained  polycrys¬ 
tals.  This  accords  with  the  premise  of  a  recent  paper  by  Gerberich  et  ah  (2003)  in 
which  size  effects  observed  in  nano-indentation  of  metallic  crystals  and  delamination 
of  thin  films  are  shown  to  scale  with  characteristic  surface-to-volume  ratios,  as  well 
as  the  increasing  normalized  strength  with  increasing  surface-area-to-volume  ratios 
observed  in  the  wire  torsion  experiments  of  Fleck  et  ah  (1994).  Potential  applications 
of  our  theory  include  stages  3  and  4  of  plastic  deformation  at  various  deformation 
rates,  as  fragmentation  into  cellular  structures  has  been  observed  in  FCC,  BCC, 
and  HCP  structures  at  large  strains  (Seefeldt,  2001),  extremely  fine-grained 
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“nanocrystalline”  materials  exhibiting  high  grain  boundary  surface  area  densities 
(Benson  et  ah,  2001;  Aleksandrov  et  ah,  2002;  Valiev  et  ah,  2002),  and  creep  condi¬ 
tions  where  subgrain  formation  frequently  occurs  (Kratochvil  and  Orlova,  1990; 
Kassner  and  Perez-Prado,  2000;  Seefeldt  and  Aifantis,  2002). 
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Appendix  A.  Kinematics  and  kinetics  of  grain  subdivision 

Here  we  illustrate  how  our  framework  may  be  applied  to  describe  evolution  of  de¬ 
fect  substructure  observed  in  ductile  (predominantly  pure  FCC)  metals  deformed  at 
low  homologous  temperatures  and  low  strain  rates  to  relatively  large  total  strains 
(Pantleon,  1996;  Hughes  et  ah,  1997,  1998,  2003;  Hughes  and  Hansen,  2000;  Hughes, 
2001).  As  shown  in  Fig.  7(a),  typically  this  substructure  consists  of  nearly  equiaxed 
“cells”  with  interiors  of  relatively  low  defect  density,  separated  by  “Incidental  Dis¬ 
location  Boundaries”  (IDBs)  across  which  lattice  misorientations  are  generally  of  the 
low-angle  variety,  for  example,  limited  to  less  than  10°  in  aluminum  and  nickel  cold- 
rolled  to  100%  effective  strain  (Hughes  et  ah,  2003).  With  increasing  applied  strain, 
the  cells  tend  to  organize  collectively  into  bands  or  “cell  blocks”,  often  elongated  in 
shape  in  a  direction  depending  upon  the  initial  texture  and  loading  conditions.  The 
cell  blocks  are  separated  by  dislocation  walls  termed  “Geometrically  Necessary 
Boundaries”  (GNBs),  across  which  misorientations  can  reach  high-angle  magnitudes 
(i.e.,  in  excess  of  15°).  Average  cell  and  cell  block  sizes  decrease  with  increasing  ap¬ 
plied  strain,  with  the  cell  block  size  generally  decreasing  at  a  significantly  faster  rate 
than  the  cell  size  (Hughes  et  ah,  1997).  Both  GNBs  and  IDBs  are  thought  to  contain 
mixed  populations  of  redundant  (i.e.,  statistically  stored)  and  non-redundant  (i.e., 


Fig.  7.  Dislocation  substructure  (a)  and  lattice  directors  (b)  in  a  deformed  single  crystal. 
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geometrically  necessary)  dislocations  (Hughes  et  al.,  2003).  It  is  also  noted  that  other 
terminology  exists  for  similar  defect  substructures;  for  example,  Kuhlmann-Wilsdorf 
and  co-workers  (Kuhlmann-Wilsdorf  et  al.,  1999)  used  the  term  “Deformation 
Bands”  to  denote  slab-shaped  volume  elements  within  which  a  reduced  number  of 
slip  systems  is  active  for  arbitrary  shape  changes  (i.e.,  fewer  than  five  systems),  with 
such  volume  elements  delineated  by  GNBs.  However,  attention  is  focused  here  on 
the  observations,  terminology,  and  hypotheses  of  Hughes,  Hansen,  and  co-workers 
detailed  in  the  aforementioned  cited  references. 

In  the  context  our  framework,  the  kinematics  of  the  above  defect  substructure  can 
be  described  as  follows.  Redundant  dislocations  and  disclinations  that  contribute  no 
net  observed  misorientation  (i.e.,  no  incompatibility)  are  grouped  collectively  into 
the  statistically  stored  dislocation  and  disclination  scalars  p's  and  fj's,  respectively. 
Geometrically  necessary  dislocations  (i.e.,  the  net  dislocation  tensor  a)  represent 
the  non-redundant  defect  component  corresponding  to  fluctuations  in  lattice  orien¬ 
tation  across  the  cell  boundaries  (i.e.,  IDBs).  The  micro-rotation  variable  Q  (in  con¬ 
junction  with  associated  net  disclination  density  tensor  6)  depicts  the  non-redundant 
dislocation  component  of  the  cell  block  boundaries  (i.e.,  GNBs).  Fractions  of  dislo¬ 
cations  and  disclinations  that  are  geometrically  necessary  may  differ,  as  they  repre¬ 
sent  different  features  of  the  microstructure:  IDBs  and  GNBs,  respectively. 

The  general  theory  developed  in  Sections  2  and  3  of  the  present  paper  places  no 
restrictions  on  the  characteristic  spacing  or  the  strength  of  the  physical  misorienta¬ 
tion  boundary  to  which  we  attribute  a  particular  type  of  defect  density.  For  example, 
we  could  alternatively  allow  dislocations  to  comprise  GNBs  and  disclinations  to 
comprise  IDBs,  or  allow  a  mixture  of  defect  types  to  comprise  each  kind  of  bound¬ 
ary.  Our  decision  here  to  model  GNBs  with  the  disclination  density  and  IDBs  with 
the  dislocation  density  stems  from  both  convenience  and  from  physical  consider¬ 
ations  in  agreement  with  models  set  forth  in  prior  literature  (Panin,  1998;  Seefeldt, 
2001).  From  the  perspective  of  calibrating  the  constitutive  model  to  physical  obser¬ 
vations,  it  is  most  convenient  to  assign  one  type  of  defect  density  (i.e.,  dislocation  or 
disclination)  to  one  type  of  boundary  (i.e.,  IDB  or  GNB)  such  that  material  param¬ 
eters  and  kinetic  relations  can  be  chosen  to  reflect  the  corresponding  physics,  for 
example  average  misorientation  angles  and  wall  spacings.  Following  Panin  (1998) 
and  Seefeldt  (2001),  we  associate  the  rotational  units  that  serve  to  accommodate 
inelastic  deformation  and  subgranular  incompatibility  seen  at  comparatively  large 
strains  with  disclinations  (GNBs),  with  geometrically  necessary  dislocations  captur¬ 
ing  the  local  misorientations  across  IDBs  arising  at  a  finer  length  scale.  From  this 
standpoint,  upon  inelastic  straining,  IDBs  (i.e.,  dislocations)  accumulate  first,  and 
then  later  organize  into  GNBs  (i.e.,  disclination  boundaries)  at  larger  applied 
deformations. 

The  micro-rotation  concept  is  illustrated  in  Fig.  7(b)  with  regard  to  lattice  vectors. 
Let  x  and  x'  be  spatial  coordinates  of  material  points  in  neighboring  cell  blocks,  with 
corresponding  lattice  director  vectors  d£,  and  (f .  Assume  that  in  elastically  unloaded 
configuration  b,  the  director  triads  are  parallel,  i.e., 

d5  =  a;.  (at) 
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Let  and  F77  denote  the  local  lattice  deformations  at  x  and  x',  respectively.  Mul¬ 
tiplying  both  sides  of  (A.l)  by  Fr'_1,  we  arrive  at  the  following  relationship  between 
d„  and  d^: 

xc) 

We  also  have,  from  (19)  and  the  assertion  of  parallel  transport  (i.e.,  V*d„  =  0), 

<  =  ((^-‘):  +  +  {F*-^dQ?aW)F%blb,  (A.3) 

where  dxc  =  x'c—xc  and  if  —  d„  w  daJ,dxb.  We  thus  see,  upon  comparing  (A.2)  and 
(A.3),  that  Q  represents  the  “error”  in  the  capability  for  the  first  spatial  gradient 
of  F7  1  to  capture  the  absolute  change  in  the  lattice  directors  over  the  small  distance 
dxc,  i.e., 


(x,c  —  xc)  )  (F^)bAb. 


(A.2) 


<  = 


(F* 


(x'c  - 


n  d  _  p&d 
z-'ca  .a 


(x'c  —  xc) 


-  (f*-1): 


(A.4) 


Equivalently,  one  can  think  of  Q  as  representative  of  the  effects  of  higher  than 
first  order  gradients  of  the  lattice  deformation  F7  on  the  spatial  variation  of  the  lat¬ 
tice  directors.  In  typical  higher-order  gradient  crystal  plasticity  theories  from  the  re¬ 
cent  literature  (cf.  Nagdhi  and  Srinivasa,  1993;  Le  and  Stumpf,  1996a;  Shizawa  and 
Zbib,  1999;  Acharya  and  Bassani,  2000;  Bammann,  2001;  Regueiro  et  ah,  2002;  Gar- 
ikipati,  2003;  Clayton  et  ah,  2004a),  only  the  first  order  gradient  of  the  elastic  or  plas¬ 
tic  deformation  (e.g.,  only  the  first  gradient  of  Fe  or  Fp  when  the  usual  decomposition 
F  =  FeFp  is  implemented)  influences  the  material  response,  through  its  manifestation 
in  the  geometrically  necessary  dislocation  density  tensor.  However,  by  including  the 
micro-rotation  variable  Q,  the  potentially  nonlinear  variation  in  lattice  directors 
(i.e.,  higher  than  first  order  gradients  of  F 7 )  between  neighboring  material  points 
is  captured  in  the  response  in  the  present  theory.  This  nonlinearity  is  manifested 
in  our  GND  density  tensor  and  in  our  disclination  density  tensor  through  Q  term 
and  its  spatial  gradient,  respectively  (see  Eqs.  (23),  (25)  and  (29)).  Physically,  this 
nonlinearity  can  be  attributed  to  the  superposition  of  lattice  rotations  across  GNBs 
between  cell  blocks  with  the  rotations  across  IDBs  between  the  cells.  Under  the 
assumption  that  the  misorientations  across  the  GNBs  remain  small,  the  disclination 
density  tensor  simplifies  to  the  linear  approximation 

4 0gd  =  Esba£dceRabce  «  28?basdced[cQe][ba]  =  2$ba^Qebafi.  (A.5) 

We  now  consider  kinetics  briefly  in  the  context  of  the  experimental  data  reported 
by  Hughes  et  al.  (1997,  1998,  2003).  Evolution  equations  applicable  in  this  context 
should  reflect  the  following  physical  phenomena: 


•  Characteristic  scaling  of  misorientation  distributions  for  IDBs  and  GNBs. 

•  Characteristic  scaling  of  size  distributions  for  cells  and  cell  blocks. 
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•  Increase  in  average  IDB  misorientation  proportional  to  (k;))  1,/2.  and  increase  in 
average  GNB  misorientation  proportional  to  (F’)2/:\  where  W  is  the  effective  mac¬ 
roscopic  plastic  strain. 

•  Reduction  in  average  cell  and  cell  block  sizes  proportional  to  ( W )~1//2  and  ( epY 2^3, 
respectively. 


Furthermore,  the  final  microstructural  configuration  attained  in  a  simulation  of 
material  behavior  should  reflect  observations  regarding  cellular  misorientation  and 
size  distributions,  and  can  be  chosen  to  reflect  a  path  of  minimum  stored  residual 
elastic  energy  (i.e.,  maximum  dissipation).  In  other  words,  for  the  same  boundary 
conditions,  more  total  dissipation  over  the  deformation  history  is  expected  for 
attainment  of  a  low-energy  configuration  consisting  of  collectively  organized  geo¬ 
metrically  necessary  and  statistically  stored  defects  than  would  be  expected  for 
attainment  of  a  plastically  homogeneously  deforming  crystal  containing  a  distribu¬ 
tion  of  redundant  (i.e.,  statistically  stored)  defects  alone.  Moreover,  greater  total  dis¬ 
sipation  is  expected  over  the  deformation  history  for  attainment  of  the  actual 
configuration  of  defects  than  for  any  other  candidate  distribution  of  these  defects. 
We  recognize  that  construction  of  specific  kinetic  relations  capturing  the  above  listed 
features  of  the  response  within  the  context  of  our  general  theory  is  a  non-trivial  task, 
one  that  will  be  undertaken  more  rigorously  in  forthcoming  works. 

Our  framework  admits  the  possibility  of  evolution  of  a  heterogeneous  deforma¬ 
tion  gradient  field  in  conjunction  with  affine  boundary  conditions  based  upon  ener¬ 
getic  arguments.  Rather  than  specify  deterministically  the  slip  rates  f  as  implied  in 
Section  3,  we  may  instead  re-write  (57)  as  incremental  relations 

di  =  dy(p‘,  p's,  fj's,  a,  0, U\  T ,  %)),  (A.6) 

where  we  have  added  dependence  upon  the  global  free  energy  of  the  entire  crystal  or 
polycrystal  (i.e.,  entire  domain  V  of  the  specific  boundary  value  problem),  with  this 
energy  denoted  by  'R  =  /  \j/dV .  Conceptually,  one  may  predict  the  evolution  of  the 
slip  increments  (A.6)  by  sequential  minimization  of  'F,  subject  to  constraints  im¬ 
posed  by  the  other  dependent  variables  and  physical  observations  (Ortiz  et  ah, 
2000).  A  non-convex  potential  i ]/  may  lead  to  lack  of  existence  of  minimum  energy 
configurations  corresponding  to  a  homogeneous  deformation  gradient  field;  instead, 
such  configurations  may  be  reached  in  the  limit  of  further  refined  microstructures 
(e.g.,  heterogeneous  deformation  gradient  fields  and  slip  rates  within  the  domain), 
as  explained  thoroughly  in  prior  literature  (cf.  Ball  and  James,  1987;  Ortiz  and  Rep- 
etto,  1999;  Carstensen  et  ah,  2002). 

To  illustrate  how  one  might  achieve  a  non-convex  energy  potential  in  the  context 
of  our  theory,  we  present  an  idealized  ID  problem.  Neglecting  for  this  simple  exam¬ 
ple  the  contributions  of  geometrically  necessary  defects,  and  assuming  fixed  temper¬ 
ature,  consider  a  non-recoverable  free  energy  i/z  depending  quadratically  upon  the 
statistically  stored  defect  densities: 


Is  1 

*  =2" 


(■ bcr)4'P2s 


rcocr 


\4~2 
)  Is 


(A.7) 
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where  cp  and  c,(  are  dimensionless  parameters.  We  let  y  denote  the  local  inelastic 
strain  in  the  crystal,  which  we  assume  evolves  monotonically,  enabling  its  use  as  a 
parameter  in  the  kinetics  equations  in  lieu  of  time.  Assume  that  the  evolution  equa¬ 
tions  (61),  mapped  to  configuration  b  and  parameterized  in  terms  of  strain,  are  of  the 
form 


dps 

0y 


A  -  Bf/S, 


djh_r 

0y 


(A.8) 


where  A,  B,  and  C  are  positive  material-dependent  constants  of  the  appropriate 
dimensions.  The  first  of  (A.8)  states  that  dislocations  accumulate  at  a  constant 
rate  minus  a  term  representing  their  trapping  and  accumulation  at  disclination 
boundaries.  The  second  simply  denotes  a  constant  rate  of  increase  in  disclination 
density  with  strain.  We  remark  that  Eqs.  (A.7)  and  (A.8),  while  perhaps  physi¬ 
cally  plausible,  present  a  drastically  simplified  rendition  of  the  general  theory 
and  are  introduced  here  for  illustrative  purposes  only.  They  are  not  necessarily 
representative  of  the  behavior  of  any  particular  material.  Assuming  initial  condi¬ 
tions  ps  =  i]s  =  0  at  y  =  0,  the  solution  of  differential  equations  (A.8)  is  found 
trivially  as 


BC  , 

Ps  =  Ay--^y ,  <is  =  cy- 


(A.9) 


Substituting  (A.9)  into  (A.7)  then  returns  the  stored  energy  as  a  function  of 
y,  i.e., 


Is  1 

*  =2" 


(bcp)4  (  Ary2  —  ABCy 3 


(rmc,^)A  C2y2 


(A.10) 


with  its  second  derivative  found  as 


02|/^ 

0y2 


( bcp )4  (  2 A2  -  6 ABCy  -  12  00  y2  j  +  2(rtnc„)4C2 


leading  to  the  conclusion  that  1 }  (y)  is  strictly  convex  when 


A1 


^'crUBCy  +  A^f. 


be, 


(ATI) 


(A. 12) 


Inequality  (  A.  12)  is  satisfied  when  C  -  0  (no  disclinations)  or  A  —  B  =  0  (no  disloca¬ 
tions),  but  is  not  satisfied  in  the  general  case  for  all  y  >  0  and  5/0.  Thus  we  see 
from  this  simple  example  that  even  though  the  stored  energy  is  convex  (quadratic 
in  this  case)  in  each  individual  defect  density  variable,  its  cross-convexity  is  not  en¬ 
sured  due  to  the  coupling  between  the  evolution  equations,  when  5/0.  Should  the 
evolution  of  local  inelastic  strain  y  be  dictated  by  energy  minimization  principles  (cf. 
Ortiz  et  al.,  2000),  the  non-convexity  of  1 }  would  be  expected  to  promote  a  spatially 
heterogeneous  distribution  of  this  strain  and  spatial  patterning  of  underlying  dislo¬ 
cation  and  disclination  densities. 
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